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ABSTRACT 

We investigate lattices of instantons and the dimension- changing transitions between 
them. Our ultimate goal is the 3D — > 4D transition, which is holographically dual to the 

phase transition between the baryonic and the quarkyonic phases of cold nuclear matter. 

O ■ 

However, in this paper (just as in [1]) we focus on lower dimensions -- the ID lattice of 

^ | instantons in a harmonic potential V oc M\x\ + M\x\ + M\x\, and the zigzag-shaped 

lattice as a first stage of the ID — > 2D transition. We prove that in the low- and moderate- 
density regimes, interactions between the instantons are dominated by two-body forces. This 
drastically simplifies finding the ground state of the instantons' orientations, so we made a 
numeric scan of the whole orientation space instead of assuming any particular ansatz. 

We find that depending on the M.2/M^/M.4 : ratios, the ground state of instanton orienta- 
tions can follow a wide variety of patterns. For the straight ID lattices, we found orientations 
periodically running over elements of a Z2, Klein, prismatic, or dihedral subgroup of the 
SU (2) /7j2, as well as irrational but link-periodic patterns. For the zigzag-shaped lattices, we 
detected 4 distinct orientation phases -- the anti-ferromagnet, another abelian phase, and 
two non-abelian phases. Allowing the zigzag amplitude to vary as a function of increasing 
instanton density, we obtained the phase diagram for the straight and zigzag-shaped lat- 
tices in the (density, M3/M4) plane. There are both first-order and second-order transitions 
between the phases, triple points, a critical point, and even a quadruple point. 

Our techniques can be applied to other types of non-abelian crystals. 
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1 Introduction 

In the ordinary N c = 3 QCD cold nuclear matter forms a quantum liquid, but for large N c 
it becomes a crystalline solid [2]. In many holographic models of QCD [3, 4] baryons are 
represented by instantons of the U(Nf) gauge theory living on the flavor branes [5], so cold 
nuclear matter corresponds to a whole crystalline lattice of instantons [6] . The geometry and 
even the dimensionality of this lattice depend on the baryon density: at medium densities, 
the instantons form a 3D lattice in the X4 = hyperplane of the holographic 4D space 
(not counting the time); this corresponds to the baryonic phase of nuclear matter. At high 
densities, the baryons spread out into the X4 dimension and form a 4D lattice; this is dual 
to the quark liquid phase of nuclear matter, or rather to the high-density quarkyonic phase 
[14] in which the quarks fill the Fermi sea, but the excitations near the Fermi surface are 



baryon-like bound states of N c quarks rather than free quarks by themselves. According to 
Rozali et al [7] , the chemical potential of the quarks in this phase is related to the thickness 
of the lattice in the X4 direction. 

Most investigations of the low-temperature holographic nuclear matter focus on the ex- 
treme high-density limit where the X4 dimension is thick and can be studied macroscopically 
(for example, see [9, 8, 7, 10]). But our main interest is in the microscopic structures of 
the instanton lattices, and especially the phase transition between the 3D and 4D lattices. 
In our previous paper [1] (with Dmitry Melnikov) we found that there is a whole sequence 
of such transitions: from a single 3D layer in the X4 = hyperplane, to two 3D layers in 
two parallel hyperplanes, to 3 layers, to 4 layers, etc., until the number of layers becomes 
too large to count individually. Alas, the 3D and 4D instanton lattices turned out to be 
too hard to analyze, so we resorted to oversimplified toy models. In the first toy model, 
we approximated the instantons as point-like charges repelling each other with 4D Coulomb 
forces; this approximation is very crude — it ignores the instantons' orientations, never mind 
the interference between the instantons — but it makes for a simple model of transitions be- 
tween lattices of different dimensions. Indeed, in this model we found that increasing density 
makes a 3D lattice of instantons in X4 = hyperplane suddenly split into two sublattices in 
the X4 = ±e hyperplanes; pictorially, this looked like the popcorn suddenly jumping up in 
the popper, so we dubbed the dimension- changing transitions the popcorn transitions. 

In our second toy model, we used actual instantons and the exact ADHM[11] construction 
for the multi-instanton system. However, we had confined the instanton centers to a 2D plane 
x\ = X2 = by making the inverse 5D gauge coupling rise steeply in the x\ and X2 directions. 
We also had the \j g\ slowly rising in the X3 direction, altogether 

-^— = N c \M(l + M 2 x\ + M' 2 (xl + x 2 2 ) + 0(M 4 Xi ;2i3 )) for M' > M. (1.1) 

Consequently, at low densities the instantons formed a ID lattice in the X4 dimension (which 
acted as the model's only flat space dimension), but for higher densities the instantons 
moved into the X3 dimension (which acted as the holographic dimension of the model). In 
that model, we saw two phase transitions in response to increasing instanton density. The 



first transition — from a straight ID chain of instantons to a zigzag-shaped chain 
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- is second-order, so the zigzag amplitude behaves as e oc \fp — p c . On both sides of this 
transition, the instantons' orientations form an anti-ferromagnetic pattern: two alternating 
orientations related by a 180° rotation. At higher densities there is another transition, but 
this time it's first-order: the zigzag amplitude e jumps discontinuously, and the instantons' 
orientations also change their pattern: the orientations of the nearest neighbors are now 
related by a rotation through angle < 110° instead of 180°. 

Presumably, at still higher densities we would have seen popcorn transitions from the 
zigzag -- which is a kind of two-layer lattice -- to a 3-layer lattice, then to 4 layers, etc., 
etc. However, calculating the net energies of such multi-layer lattices turned out to be too 
hard, so we stopped at the zigzag. 

Technically, the main difficulty in working with more complicated instanton lattices is 
setting up the ADHM construction [11] for the infinite number of instantons. Indeed, the 
ADHM construction for A instantons involves 4 non-commuting A x A matrices T^ n whose 
off-diagonal elements follow from non-trivial constraint equations; for A — > 00 solving these 
equations becomes very hard. For a simple ID lattice of instantons whose orientations 
are related by commuting SU(2) symmetries, the ADHM constraints were solved by Kraan 
and Baal [12], and we have re-derived and used their explicit formula for the instanton 
number density profile I(x) in [1]. For a zigzag-shaped instanton chain with link-periodic 
instanton orientations, we found an exact solution of the ADHM constraints in terms of 
Fourier transforms of some rather messy functions, but we could not use it to derive an 
exact formula for the I(x); instead, we had to use the small-instanton-size approximation. 
And it gets worse for more complicated instanton lattices and/or orientation patterns: For 
the simple 2D and 3D lattices and purely abelian orientation patters the ADHM constraints 
have known formal solutions 1 - - which alas are way too complicated for any practical use 

- while for the non-abelian orientation patterns there are no known solutions at all. This is 
particularly unfortunate since the lowest-energy 3D configuration for low densities is likely 
to be an FCC lattice with non-commuting orientations of the instantons [13] . 



*A Nahm transform can map the ADHM constraints to an electrostatic- like problem in 2D or 3D. 



In this paper we follow a shortcut around the ADHM construction -- and also around 
the more difficult part of the energy calculation. We assume from the beginning that the 
instanton density is not too high, so the distances between instantons are much larger than 
the instantons' radii. In §3 we show that for all such multi-instanton systems, interac- 
tions between the instantons are dominated by the two-body forces. The irreducibly-three- 
body forces, 4-body forces, etc. also exist, but they are suppressed by powers of the small 
(radius/distance) 2 ratio compared to the two-body forces. Consequently, we obtain a man- 
ageably simple approximate formula for the net energy of a multi-instanton system in terms 
of the instantons' positions and orientations. 

Thanks to this formula, we do not need to assume (as we did in [1]) that in a lattice of 
instantons, their relative orientations have the same symmetries as their relative positions. 
Instead, we may use numerical methods to minimize the net energy with respect to all the 
orientations treated as independent parameters, without any symmetry assumptions. The 
algorithm we use in this paper is fairly simple: We set up a large lattice (200 sites is more 
than enough for 1 dimension) of fixed geometry, with a free SU(2) matrix y n at each site 
encoding the n th instanton's orientation (for Nf = 2). We start with a completely random 
set of the initial y n and make them evolve along the steepest descent of the net energy 
function until they converge to a local minimum. Then we repeat the process many times 
with different random sets of initial y n . This gives us a short list of the local minima, from 
which we select the deepest -- which is presumably the global minimum of the net energy 
function. This algorithm has convergence problems near the phase boundaries — where two 
or more local minima are almost degenerate -- but it is good at finding all the orientation 
patterns that might be ground states for some combinations of background parameters and 
lattice geometries. And once we know all these patterns, we fit them all with common ansatz, 
calculate the net energy analytically as a function of the ansatz's parameters, and then we 
minimize it with a much higher precision to get the accurate phase diagram of the instanton 
lattice. 

In this paper we put the instantons in a harmonic potential that rises at different rates 
along three out of four space coordinates; in terms of the 5D gauge coupling, 
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gf(^l,X 2 ,X 3 ,X4) 



N C XM x (l + Mlx\ + M\x\ + M\x\ + 0(M 4 x^ 4 )) (1.2) 



with 3 independent parameters M 2 ,M 3 ,M^. Here we use different coordinate axes than in 
[1]: the ID straight chain — or the long axis of the zigzag — runs along the X\ axis instead of 



the X4. But the important difference is allowing for M 3 7^ M 4 and M 2 ~ M 3 ~ M 4 instead of 
M 3 = M 4 = M' ^> M = M 2 ; it is this difference which leads to the wide variety of instanton 
orientation patterns we shall see in this paper. 

Specifically, we shall see that for the ID instanton lattice in the almost-holographic 
background M 2 ,M 3 <C M4 (where M 2 and M 3 help to line up the instanton centers along 
the x\ axis but do not affect the two-body forces between the instantons), there is a large 
degenerate family of ground states for the instanton orientations y n G SU(2)/Z 2 (for Nf = 2). 
This family includes many periodic patterns in which all the y n belong to a finite subgroup of 
the SU(2)/Z 2 : the anti-ferromagnetic pattern in which y ddn and y cvcnn belong to 2 elements 
of a Z 2 subgroup, the period = 4 pattern in which the y n span the Klein group Z 2 x Z 2 , and 
two series of period = Ik patterns in which the orientations span a prismatic group Z& x Z 2 
or a dihedral group D 2 k- There are also many link-periodic patterns in which the relative 
orientations y n y n+ \ of neighboring instantons are periodic but the orientations y n themselves 
are not periodic. But the vast majority of the degenerate ground states are not periodic at 
all. 

Increasing the M 2 and M3 parameters lifts the degeneracy. For M 2 <C M3 ~ M4, the 
ground state of the ID lattice is anti-ferromagnetic, while for M 2 ~ M 3 ~ M4 the ground 
state is periodic with a longer period or link-periodic, depending on the M 2 /M% ratio. Specif- 
ically, for M 2 = M3 the instanton's orientation in the ground state span the Klein group, 
for other rational M\jM\ they span a dihedral group, while the irrational M^/M^ favor the 
link-periodic patterns. 

For the zigzag-shaped chains of instantons we need M 2 <^ M 3 ~ M 4 to make sure the 
transition from a straight chain to a zigzag happens when lattice spacing ^> instanton radius 
- otherwise the two-body-force approximation would not work - - but we may vary the 
M3/M4. In this setting, the ground state of instanton orientations may follow one of four 
different patterns, depending on the zigzag amplitude and on the M3/M4 ratio. Besides the 
anti- ferromagnetic and the abelian <p < tt patterns seen in [1], there are two non-abelian 
link-periodic pattern in which the 2/^2/n.+i "twists" for even and odd links of the lattice do 
not commute with each other. The non-abelian patterns appear only for M 3 < M 4 , that's 
why we did not see them in [1]. 

Focusing on the "popcorn transition" from a straight chain of instantons to a zigzag- 
shaped chain (which is a two-layer lattice), we shall see that the thermodynamic order of 
this transition depends on the M3/M4 ratio: the transition is first-order for M3 < O.63M4 
but second order for M3 > O.63M4. Likewise, when we further increase the instanton density, 
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the number and the orders of the subsequent phase transitions between different orientation 
phases of the zigzag-shaped chain also depend on the M 3 /M^ ratio. 

This paper is organized as follows: For the reader's convenience, in §2 we review the 
basic features of holographic nuclear matter: effects of the large N c and large A limits, the 
generalized Sakai-Sugimoto model of holographic QCD, realization of holographic baryons 
as instantons of the flavor gauge symmetry, basic features of the holographic nuclear forces, 
and our basic approach to the multi-baryon systems and the bulk nuclear matter. In §3, 
we argue that the interactions between multiple instantons include 2-body forces as well as 
irreducibly-3-body force, 4-body forces, etc., but when the distances between the baryons 
are much larger that the instantons' radii, the two-body forces dominate the net energy, 
while the 3-body,. . . , forces are suppressed by powers of (radius/distance) 2 <C 1. First we 
give a heuristic argument, and then we back it up with a technical proof for Nf = 2; we 
also calculate the precise form of the net two-body force and its dependence on the two 
instantons' orientations. 

In §4 we use the two-body-force approximation to study ID lattices of instantons. We 
start with the M 2 ,M 3 -C M4 background and show that there is a very large degenerate 
family of ground states of the instantons' orientations. Then we move on to the M 2 ~ M 3 ~ 
M4 backgrounds; this changes the precise form of the two-body force in a way that lifts the 
degeneracy between the orientation states. Consequently, the ground state of the lattice has 
a specific periodic or link-periodic pattern of orientations, depending on the M 2 /M 3 ratio. 

In §5 we study zigzag-shaped instanton chains in the M 2 <C M 3 ~ M 4 backgrounds. After 
a few preliminaries, we use numerical methods to find the lowest-energy state of instanton 
orientations without presuming any particular symmetries. We repeat the calculation for 
50 different lattice geometries (e/D = 0.01,... 0.99 by 0.02) and 50 values of the M3/M4 
ratio; the results are presented on the phase diagram 5 on page 63. Our numeric code 
is not very accurate and has convergence problems near the phase boundaries, but it is 
good at identifying the symmetries of the lowest-energy patterns. Using these symmetries, 
we construct a 3-parameter ansatz that includes all the preferred patterns, then calculate 
the net energy as analytic function of the parameters and minimize it. For the reader's 
convenience, the energy calculation and minimization is described in the Appendix, while 
the resulting phase diagram is presented in figure 6 (page 67) in the main body of §5. 

In the last part of §5 we treat the zigzag amplitude as a dynamic parameter, vary the 
instanton density, and study the popcorn transitions from the straight chain to the zigzag 
and between different orientations phases of the zigzag. We show that the net number of 
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such transitions -- 1, 2, or 3 -- and their orders depend on the M3/M4 ratio. The overall 
phase diagram in the density v. M3/M4 plane is presented in figure 7 on page 71. 

Finally, in section 6 we summarize our results and give a short list of open questions we 
hope to address in the near future. 

2 Review of Holographic Nuclear Physics 

In this section we review the basics of holographic nuclear physics that we shall need later 
in this paper. The holography here means gauge-gravity duality, and the semiclassical de- 
scription of the gravity side of the duality requires the limits of large iV c and large 't Hooft 
coupling A = #ym x N c [15]. In the first two subsections §2.1-2 we discuss the general con- 
sequences of these limits for the nuclear physics. In §2.3 we get more specific and review 
the Sakai-Sugimoto model of holographic QCD. There are many other models which differ 
in countless details - - and none of them is exactly dual to the QCD, - - but their basic 
features are rather similar, and the Sakai-Sugimoto model serves as a prototype. In §2.4 we 
introduce the baryons and explain their holographic duality to the instantons of the flavor 
gauge symmetry. We also review the interactions and the self-interactions of such instantons 
in the large A limit and explain why they have small radii a oc A -1 / 2 . Finally, in §2.5 we 
address the multi-instanton systems dual to nuclear matter and explain how the interactions 
between the baryons appear in at the first order in the 1/A expansion. We also explain the 
technical difficulties in handling infinite crystals of baryons and preview a shortcut around 
those difficulties that we shall take in §3. 

2.1 Nuclear matter in the large N c limit 

QCD perturbation theory in the large N c limit is dominated by planar gluon diagrams while 
contributions of the non-planar diagrams and of the quark loops are suppressed. At zero 
temperature and chemical potential below the mass of the baryon, the large N c QCD is a 
theory of free mesons and glueballs whose interactions are suppressed by powers of N c (except 
for interactions with the baryons). The masses of mesons and glueballs are 0(Aqcd)> and 
they scale as N®. On the other hand, the baryons - - which are made out of iV c quarks 
- have masses of order N c x Aqcd, so their relative abundance in thermal equilibrium is 
exponentially suppressed. The interaction energy between baryons also scales as N c . 



The T — li phase diagram of the large-iV c QCD is believed to look like figure 1 below. 
(For a review see [14] and references therein, also [22].) We assume the 't Hooft's iV c — > oo 
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Figure 1: Conjectured phase diagram for QCD at large N c [14] 

limit in which the number of flavors remains finite. Consequently, the dynamics of the 
theory is dominated by the gluons, the quarks are sensitive to the gluonic background, but 
the backreaction from the quarks to the gluons is suppressed by Nf/N c <C 1. Thus, at 
lower temperatures there is confinement but for increasing temperature there is a first order 
transition to the deconfined phase. The transition temperature Td is at the Aqcd scale and 
it's almost independent of the quark chemical potential ii q = fXb/N c (as long as /x q does not 
grow with the N c ). 

It is not clear whether for N c — > oo the deconfining phase transition coincides with 
the chiral symmetry restoration for the light quarks. Several field theory arguments -- for 
example [16, 14] - - suggest that for fx q = the two transitions should happen at the same 
point. However, these arguments do not work for fx q > [14], and there are other arguments 
for the existence of confined but chirally restored phases and/or deconfined phases where 
the chiral symmetry remains broken. Moreover, some holographic models — for example [4] 
- have deconfined but chirally broken phases even at zero chemical potentials. 

For temperatures below the deconfining transition td and baryon chemical potentials 
below the baryon mass - - or equivalently for ii q < m q = Mb/N c - - the thermal state of 
the theory is a gas of glueballs and mesons with almost no baryons or antibaryons. But at 
fx q m m q there is an abrupt phase transition to the bulk nuclear matter with finite baryon 
density. But unlike the ordinary nuclear matter which is in the quantum liquid state, the 
large N c nuclear matter is crystalline solid since the ratio of kinetic energy to potential energy 
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decreases with N c . Indeed, the potential energy of baryon-baryon forces scales like JV C ; more 
precisely [17], in the large N c limit the two-baryon potential becomes 

V ~ N c x A c (r) + N c x A s (r) (h • I 2 ) (J x • J 2 ) + 

+ N c x A T (r) (Ii • I 2 ) (3 (n • 3 1 ) (n • J 2 ) - (J x • J 2 )) + O (1/N C ) , (2.1) 

for some ^-independent profiles Ac(r), As(r), and Ar{r) for the central, spin-spin, and 
tensor forces; their overall magnitudes are A ~ Aqcd for r ~ 1/Aqcd- Classically, this 
potential tries to organize the baryons into some kind of a crystal where the distances between 
neighboring baryons do not depend on the N c while the binding energy (per baryon) scales 
like N c x O(Aqcd)- in quantum mechanics, the baryons in such a crystal behave like atoms 
in ordinary crystals — they oscillate in their potential wells with zero-point kinetic energies 

K *9£dL (2.2) 

2m B d 2 N c d 2 ' V } 

where d ~ 1/Aqcd is the N c - independent diameter of the potential well. Therefore, at zero 
temperature the ratio of kinetic energy to the potential energy scales like 

(2.3) 

and becomes very small for large N c . At higher temperatures the kinetic energies of baryons 
become larger, K ~ T, but in the confined phase we are limited to T < Td ~ Aqcd and 
hence K < O(Aqcd)- Consequently, the kinetic to potential energy ratio scales like 

(2.4) 

- which is larger than (2.3) but still becomes small in the large A^ c limit. Consequently, 
for large A^ c neither zero-point quantum motion nor thermal motion of baryons can melt the 
baryon crystal, so nuclear matter remains solid all the way to the deconfining temperature. 

Before holography, the best models for the large-Af c nuclear crystals were lattices of 
skyrmions. In this framework, Igor Klebanov [18] had found a curious phase transition from 
a lattice of whole skyrmions at low chemical potentials to a denser lattice of half-skyrmions 
at higher potentials. (According to [19], the whole- skyrmion lattice at low potentials has 
FCC geometry while the half-skyrmion lattice at higher potentials is simple cubic.) In the 
half-skyrmion-lattice phase, the order parameter for the chiral symmetry breaking vanishes 
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after space averaging, so in QCD terms the transition is interpreted as chiral symmetry 
restoration at high /i q . 

In QCD with N c = 3 and two massless flavors there is a similar chirally-symmetric phase 
at low T and high fi q . This phase is a quark liquid rather than a baryon liquid (the quarks are 
no longer confined to individual baryons) and there is a condensate of quark pairs making 
this liquid a color superconductor. But for large N c the situation is more complicated: 
there is no color superconductivity, and there is no deconfinement for T < Tj. Instead, the 
dense cold nuclear matter forms a phase which combines the features of the baryonic and 
quark phases: the quarks fill up a Fermi sea, but the interactions near the Fermi surface are 
strong, so the excitations are not free quarks or holes but rather meson-like quark-hole pairs 
or baryon-like states of N c quarks. MacLerran et al have dubbed this phase "quarkyonic" . 

For fi q ^> Aqcd, the interior of the Fermi see is chirally symmetric, but near the Fermi 
surface the symmetry is broken by the chiral density waves [20]. (Although Son and Shuster 
[28] argue that such waves develop only for very large N c > 10 3 .) To be precise, the chiral 
density waves mix the chiral symmetry of the quarkyonic phase with the translational sym- 
metry rather than simply break it. Averaging over space restores the chiral symmetry, just 
like it happens for the lattice of half-skyrmions. 

On the other hand, for \i q just above m q = Mb/N c , the baryonic crystal has a completely 

(c) 

broken chiral symmetry. Thus, at some critical /j, q = m q x a few, there should be a chiral 
symmetry restoring phase transition from the baryonic crystal to a distinct quarkyonic phase. 
In holography, this transition is believed to be dual to the "popcorn transition" from a 3D 
instanton lattice to a 4D lattice. 

2.2 Effect of the large A limit on the holographic nuclear matter 

In holography, the semiclassical description of the gravity side of the gauge-gravity duality 
in terms of metric, fluxes, and branes requires the limits of large N c and also large 't Hooft 
coupling A = NcQym- I n the large A limit, the baryons become very heavy: in units of the 
mesonic mass scale M ~ Aqcd, the baryon mass is M& ~ XN C M. However, the interactions 
between the baryons do not grow with A: even for two baryons right on top of each other, 
the repulsive potential between them is only V ~ N C M ~ M^/X. At larger distances, the 
forces are even weaker since the hard-core radius of a holographic baryon shrinks with A as 
R b ~ M~ 1 A -1 / 2 . Outside this radius, the repulsive potential decreases as 1/r 2 until r ~ M -1 , 
at which point it becomes dependent on the meson mass spectrum of a specific holographic 
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model: In some models, the potential becomes attractive for r > M l while in others it 
remains repulsive at all distances. The overall picture is shown in figure 2: 



0{M b /\) 



Model 




A 2 



OiX-^/M) 



0{l/M) 
Figure 2: Two-body nuclear potential in holographic QCD. 



Since the nuclear forces are so weak in the holographic QCD, all transitions between 
different phases of the cold nuclear matter happen at chemical potentials /if, very close to the 
baryon mass: just below Mb we have glueball/meson gas (or vacuum for T = 0), while just 
above Mb we have dense quark matter. To see the baryonic matter phase (or any other inter- 
mediate phases) we need to zoom into the //& ~ M b region. Figure 3 (on the next page) illus- 
trates this point: to see the baryonic-matter phase between the vacuum and the quark-matter 
phases on the plot of baryon density p as a function of the chemical potential /i, we need to 
zoom into the narrow range \x — Mb = 0(Mb/X). The figure also shows that the thermody- 
namic order of the phase transition between the vacuum (or the meson/glueball gas for T > 
0) and the baryonic matter depends on the sign of the long-distances nuclear force. If the force 
becomes attractive at long distances, then bulk baryonic matter exists at zero external pres- 
sure and has /i = Mb — binding energy. Consequently, the transition from the vacuum (or gas) 
to the nuclear matter (in the form of a baryonic crystal) is first-order as shown on figure 3(a). 
On the other hand, if the nuclear forces are repulsive at all distances, then the bulk nuclear 
matter does not exists except at positive external pressures, and its chemical potential must 
be \i > Mf,. Moreover, \x raises monotonically with the pressure and the density, so the tran- 
sition from the vacuum to the bulk nuclear matter is second-order as shown on figure 3(b). 
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(a) 



(b) 




Figure 3: Density as a function of the chemical potential in the large A regime for attract- 
ing (a) and repelling (b) baryons. The nuclear matter phase is confined to a narrow window 
of the order Afi ~ Bg/X. In the naive diagram the transition occurs directly from the 
no-baryon to a quark phase. 

In this article we do not wish to focus on a particular model of holographic QCD, so we 
make no assumptions about the long-distance nuclear forces. Consequently, we cannot say 
anything specific about the transition from the vacuum to the baryonic phase. Instead, we 
focus on the transition from the baryonic phase to the quarkyonic phase, which correspond 
in the holographic picture to changing the lattice geometry of the instanton crystal, from 
a 3D lattice, through a sequence of intermediate steps, to a 4D lattice. Or rather, that's 
the goal of our program; in this article we shall focus on a simplified problem, namely the 
transitions between ID and 2D instanton lattices. 



2.3 Sakai— Sugimoto model as a prototype of holographic QCD 

Some gauge theories — especially the M = 4 SYM -- have exact holographic duals, where 
both sides of the duality follow as IR limits of the same string-theoretical construction, and 
all the extra degrees of freedom are superheavy. Alas, the Quantum Chromo Dynamics — or 
even the large N c limit of QCD — is not so lucky: it either does not have an exact holographic 
limit, or we have not found it yet. Instead, there is a large number of "holographic QCD" 
models, which are dual not to the QCD as such but rather to some QCD-like theories with 
a lot of extra junk that the real QCD does not have. Such models make for "spherical horse 
in a vacuum" approximations — good enough to get qualitative understanding, but not too 
accurate to make useful numerical predictions. 
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In a moment, we shall review the Sakai-Sugimoto model [3] of holographic QCD. That 
model should not be taken too seriously — its predictions for the hadronic spectrum are at 
best so-so, and for the nuclear forces are even worse, - - but it serves are a type specimen 
for a large class of models. Here are some general features of this class: 

• The construction starts with N c coincident D-branes, which span the Minkowski space x 
a compact cycle. Typically, one uses D3 branes at a singularity, or D4-branes wrapping 
a circle, although some constructions use D5 branes on a 2-cycle, etc. For weak 't Hooft 
couplings A = N c x g Y M ~ N c x g str <C 1, the open strings between the branes give rise 
to the gluons of the U(N C ) gauge theory (plus a lot of extra junk). 

• For A ^> 1, the D-branes merge into a big fat black brane. All we see outside the 
horizon is a warped space-time geometry and the Ramond-Ramond fluxes induced 
by the conserved charges of the D-branes. Modes of the bulk-field excitations in this 
geometry are dual to the QCD glueballs. Or rather, some modes are dual to the 
glueballs while other models are dual to the extra junk that the real QCD does not 
have l^. 

• To make the hadrons (mesons and baryons), add Nf flavor D-branes, usually D5, D6, 
D7, or D8. At weak A, the open strings connecting the color branes to the flavor branes 
give rise to the quarks and the anti-quarks. 

• In the holographic limit N c — > oo and A — > oo, but the Nf stays finite and the Nf x ^ str 
remains weak. Consequently, the flavor branes remains D-branes (rather than merge 
into a black brane of their own) , but now they live in a warped geometry with fluxes 
induced by the N c color branes. 

• The open strings between the flavor branes give rise to Nj vector and scalar fields 
living on those branes. The 4D modes of these vector and scalar fields are dual to the 
QCD mesons. 

• The YM instantons of the vector fields on the flavor branes are dual to the QCD 
baryons, see §2.4 for more details. 

And now let us review the Sakai-Sugimoto model [3] in a little more detail. The color 
sector of the model is the Witten's model [23]: N c D4 branes spanning the Minkowski 
space x a circle; the circle has antiperiodic boundary conditions for the fermions, which breaks 
the Af = 4 SUSY down to M = 0*. In the field theory limit A < 1 => A QCD < M KK 
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the effective low-energy theory is pure U(N C ) Yang — Mills, but in the holographic limit 
A ^> 1 everything happens right at the Kaluza-Klein scale Mrk = 1 /circle radius, so the 
YM glueballs end up with similar O(Mkk) masses to a lot of non-YM junk. 

On the gravity side of the duality, the D4 branes merge into a black brane which warps the 
10D metric. Instead of flat 10D spacetime, we now have a warped product of R 3,1 Minkowski 
space, the S 4 sphere (originally surrounding the D4 branes), and a 2D cigar spanning the 
radial direction (_L to the branes) and the S 1 circle. The radial coordinate u runs from 
ma > to infinity; at u\ the S* 1 circle shrinks to a point, hence the cigar. Altogether, we 
have warped metric, the 4-form flux, and the dilaton's gradient according to 



ds 2 



\ 3/ 2 

7/ \ 



RdaJ 



-dt + bijdx % dx 3 + f(u)dx 



+ M" 


- du 2 
TTT + u2dn l 


r J "Y" 





F 4 = 3irt s N c x volume_form(S 4 ) , e* = g s [ —— , (2.5) 

\Rd4J 

where x 4 was the coordinate along the S 1 circle and now is the polar angle on the cigar, 

R 3 m = 7rg s elN c , /(«) = 1 - (^) 3 , (2.6) 

£ s = va 7 is the string length scale, and g s is the string coupling. The U\ -- the minimal 
value of the radial coordinate u at the tip of the cigar — is related to the original radius R 
of the S 1 circle as 

2lfl ^(^y /2 . (2 . 7) 

3 V M A / 
The same radius R also controls the 4D Yang-Mills coupling and hence the 't Hooft's coupling 
A. Analytically continuing from A <^ 1 to A ^> 1, we have 

A = gt D x N c = || x N c = ^ x N c . (2.8) 

To add the flavor degrees of freedom to the model, Sakai and Sugimoto added Nf D8 
branes and Nf anti-D8 branes. On the field-theory side of the holographic duality (i. e., for 
A <C 1), they span all space coordinates except the x 4 coordinate along the S 1 circle. At the 
intersections of D8 and D4 branes, the open strings give rise to the massless quarks; likewise, 
at the intersections of the anti-D8 branes and the D4 branes we get massless antiquarks. 

On the holographic side A ^> 1 of the duality, the exact solution for the flavor branes 
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interacting with the warped metric and fluxes is not known, but for Nf <C N c and Nf x 
g s C 1 we may use the probe approximation: The flavor branes seek the lowest-action 
configuration in the background metric (2.5), while their back-reaction upon the metric is 
neglected. Consequently, at low temperatures (below the deconfinement transition) 2 , the 
flavor branes span the Minkowski space xS 4 xa line on the cigar; the exact shape of this 
line follows from minimizing the branes' action, but its topology follows from the cigar itself: 
Since the D8 and the anti-D8 branes cannot reach all the way to the origin u = 0, they 
must reconnect to each other and form U-shaped configurations as shown on figure 4. The 





Figure 4: The figure on the right is the generalized non-antipodal configuration. The figure 
on the right describes the limiting antipodal case L 
u = u A . 



irR, where the branes connect at 



reconnection is the geometric realization of the spontaneous chiral symmetry breaking: the 
separate stacks of Nf D8 and Nf anti-D8 probe branes give rise to the U(Nf) L x U(Nf) R 
gauge symmetry, which corresponds to the U{Nf)i x U(Nf) R global chiral symmetry in 
4D. But when the D8 and anti-D8 branes reconnect, there is only a single stack of Nf U- 
shaped branes and hence only one unbroken U(Nf) symmetry. Thus, the chiral symmetry 
is spontaneously broken, 

U(Nf) L xU(Nf) R ^U(N f ). (2.9) 

As shown on figure 4, the U-shaped profiles of reconnected branes depend on one param- 
eter — the asymptotic distance L between the D8 and the anti-D8 branes along the S 1 circle 
for u — > oo. For L = nR the branes form the antipodal configuration in which the branes 
remain at opposite points on the circle for all u, all the way from u — > oo down to « = % 



2 At higher temperatures — above the deconfining transition — the background metric has different 
topology, and the flavor branes also have different shapes, see Aharony et al [4] for details. 
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where the branes reconnect; this is the original configuration of Sakai and Sugimoto. In the 
more general version of the model [4] we allow for the L < nR non-antipodal configurations. 
In such configurations, the distance between the branes in the X4 direction depends on u — it 
becomes smaller for smaller u — and eventually the branes reconnect at uq before they reach 
the bottom of the cigar. The uq/u^ ratio may be used to parameterize the non-antipodal 
configurations instead of the L/R. 

The Uo/ux or L/R parameter of the Sakai-Sugimoto model does not correspond to any 
adjustable parameters of the real-life QCD. Unfortunately, this parameter affects many phys- 
ical properties of the model. For example, for {L/R) > 0.97 the deconfinement and the 
restoration of chiral symmetry happen at the same temperature, but for (L/R) < 0.97 they 
happen at different temperatures and the model has an intermediate deconfined but chirally 
broken phase [4]. Also, in the antipodal model the central nuclear forces are purely repulsive, 
while the non-antipodal models give rise to both repulsive and attractive nuclear forces [2] 
(but unfortunately the net force remains repulsive at all distances). 

The low-energy dynamics of the flavor degrees of freedom living on the D8 branes is 
governed by the effective action comprising the Dirac-Born-Infeld (DBI) and Chern-Simons 
(CS) terms, 

S = S DBI + S cs ■ (2.10) 

The DBI action is 

Sdbi = T 8 [ d 9 x e-* Str ( y/- det(g mn + 2Tia'J mn )^ , (2.11) 

D8+D8 

where T$ = (2 , k)~' p £~ p ~ 1 is the D8-brane tension, g mn is the nine-dimensional induced metric 
on the branes, T mn is the U(Nf) gauge field strength, and Str denotes the symmetrized trace 
over the flavor indices. In the limit of fixed brane geometry and weak gauge fields, the DBI 
action reduces to Yang-Mills, 

SdbiIF] = const + S YM [T] + 0(J^ 4 ). (2.12) 

Furthermore, the low-energy field modes we are interested in are constant along the S A sphere 
and the vector fields' directions are _L S 4 . Consequently, we are going to dimensionally 
reduce the flavor gauge theory down to 5 dimensions: the 4 Minkowski dimensions re 0,1 ' 2,3 , 
plus one coordinate z along the U-shaped line on the cigar. We find it convenient to choose 
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a particular z coordinate that makes the 5D metric conformal 



A( z ) x (-dt 2 + rfx 2 + dz 2 ), A{z) = (^L) . 



ds 2 = A(z) x (-dt 2 + dx 2 + dz 2 ), A(z) = -^- ; (2.13) 



although its relations to the u and x 4 coordinates of the cigar follow from the differential 
equations 



du(z) = u*(z)f(u(z)) - u 8 J(u Q ) dx A (z) = U^y/fjup) 

dz ~ y u 5 (z)R 3 D4 ' dz " u\z)f{u{z)) ' l ' j 

In the (x°,x l ,x 2 ,x 3 , z) coordinates, the 5D YM action for the flavor gauge fields becomes 

Sym ~ \d A x dz 2 tr(-^AfAr) (2-15) 

for 

1 = ^ M KK «(£) , . 

2<? YM (z) 2167T 3 MA • l • ^ 

Near the bottom of the U-shaped flavor branes — which is the only region we are going to 
care about in this article — 



1 = Ay_MKK / 8C^5 + \ where «, 

2 ^ym(^) 216tH V 9 C / m a 

(2.17) 

The Chern-Simons term arises from the couplings of the gauge fields on the D8 brane to 
the bulk Ramond-Ramond field. In 9 dimensions 

S cs =T 8 x I CgAtre 2 ™'^, where F 4 = dC 3 . (2.18) 

D8+D8 

After integrating over the S A and dimensionally reducing to 5D, the Chern-Simons term 
becomes 

Scs = 2& L tr (AF2 ~ M ' 3 ^ ~ ^ 5) • (2 - i9) 

In a particularly interesting case of 2 flavors, it is convenient to separate the U(2) gauge 
fields Am into their SU(2) components Am and the U(l) components Am- In terms of these 



components, the Chern-Simons action becomes 

*» = I^/^tr^ + ^JaaF>. (2,0) 

We shall see in a moment that the baryons and the multi-baryon systems have strong self- 
dual £77(2) magnetic fields F^. 3 Thanks to the first term in this Chern-Simons action, the 
instanton number density 

I(x,z) = 32 g ^ (2.21) 

acts as electric charge density for the abelian field A ; the net electric charge of an instanton 
is Q e i = N c /2. 

Besides the U(Nf) gauge fields, the effective low-energy 5D theory also contains the scalar 
fields <3> a (x, z) describing the small fluctuations of the D8 branes in the transverse direction. 
For Nf branes, the scalars form the adjoint multiplet of the U(Nf) gauge symmetry. The 
action for the scalar fields follow from the DBI action for the embedded metric g mn of the 
fluctuating branes. For the <&(x, z) fields normalized to have similar kinetic energies to the 
vector fields, the scalar action looks like 

Scalar = f#X f dz —^—- tl((D M <$>f + V (<$>) 

J J 2 9ym{ z ) v 

+ ^ fd 4 xfdzC(z) x ti(^ MN J MN ) + ■■■ . (2.22) 

The details of the scalar potential V(<&) = m 2 (z) x $ 2 + a(z) x $ 4 + • • • need not concern 
us here, what's important is the second term describing the backreaction of the gauge fields 
on the brane geometry. In the antipodal Sakai-Sugimoto model C(x) = and there is no 
backreaction because of a geometric symmetry, but in the non-antipodal models C(z) ^ 
and the scalar fields $ induced by the vector fields of the baryons lead to attractive nuclear 
forces [2]. The ratio of these attractive forces to the repulsive forces mediated by the abelian 
electric fields depends of the C(z) profile of the interaction term (2.22). For the Sakai- 
Sugimoto models 

attractive ~. . 1 — C~ 3 / Un \ 1 , 

C 2 (z) = -^— x — \ < - < 1, (2.23) 



repulsive 9 \u{z) J 9 



3 In our notations, the spacetime indices 0,1,2,3, z of the effective 5D theory are labeled M,N,... = 
0, 1, 2, 3, z while the space indices 1, 2, 3, z of the same theory are labeled /a, v, . . .. When we need the 9D 
indices for the whole D8 brane, we use m,n, . . .. 
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so the net force is repulsive. 

To see how that works, let us focus on the baryons in the Sakai-Sugimoto and other 
models of the holographic QCD. 

2.4 Baryons as flavor instantons 

In the old hadronic string, the baryons were made out of Y-shaped configurations of 3 open 
strings connected to each other at one end; the other end of each string is connected to a 
quark. To realize this picture in holographic QCD, we need a baryon vertex (BV) -- some 
object connected to N c open strings. The other ends of the strings are connected to the flavor 
branes and act like the quarks; this gives the baryon its flavor quantum numbers. Witten 
had constructed the baryon vertex for the AdSs x S* 5 model from a D5 brane wrapping the 
S 5 sphere [24]; the generalized versions of this construction [25] use D p branes wrapping 
compact cycles carrying 0(N C ) Ramond-Ramond fluxes. 

In the Sakai-Sugimoto version of this construction, the baryon vertex is realized as a D4 
brane wrapped around the S 4 sphere (but localized in all other dimensions except the time). 
The S 4 carries N c units of the F^ Ramond-Ramond flux, 



so the Chern-Simons coupling of this flux to the U(l) gauge field B living on the D4 brane 
acts a N c units of the net electric charge for the B : 

T 4 f C 3 A e 2wa ' dB = N c f B dt . (2.25) 

</D4 ./time 

In a compact space like S 4 , the net electric charge must vanish. To cancel the charge (2.25) 
we need to connect the D4 brane to open strings. The back end of an oriented open string 
has electric charge — 1, so we must connect the D4 brane with N c such strings; their front 
ends connect to the D8 flavor branes (since the strings do not have any other place to end) 
and act as iV c quarks. 

We may put the D4 brane anywhere in space and anywhere on the cigar. However, the 
S 4 volume increases with the u coordinate, so the lowest-energy location of the D4 is the 
cigar's tip u = U\. At other location, the brane feels a gravity-like force pulling it down to 
the tip. However, the strings connected to the baryon vertex pull it towards the flavor D8 
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branes; in the non-antipodal models the D8 branes do not reach the cigar's tip, so the strings 
pull the baryonic vertex up from the tip towards the lowest point uq of the flavor branes. 
The 'tag of war' between the upward and downward forces on the baryon vertex determines 
its ultimate location. In some models, the forces reach equilibrium for the baryon vertex 
hanging on strings below the flavor branes [35], while in many other models — including the 
non-antipodal Sakai-Sugimoto -- the string forces win and pull the baryon vertex all the 
way to the lowest point Uq of the flavor branes [30]. 

In all such models, the baryonic vertex is a D p brane completely embedded in a stack 
of D p+4 flavor branes, so it is equivalent to zero-radius Yang-Mills instanton of the U{Nf) 
gauge symmetry on the flavor branes, and for Nf > 1 it may be smoothly inflated to a 
finite-radius instanton [31]. In p + 5 dimensions of the flavor branes, this instanton is a fat 
Dp-brane wrapping some compact cycle, but once we dimensionally reduce to 5 dimensions, 
the instanton becomes a finite-size particle. Thus, in the low-energy effective 5D theory of the 
holographic QCD, a baryon is realized as a finite-size instanton of the U(Nf) gauge theory. 

Our research program of studying the instanton crystals as holographic duals of the 
cold nuclear matter stems from this realization of baryons as instantons in 5D, but it does 
not depend on any specific details of the Sakai-Sugimoto model. We use that model as 
a prototypical example, but any other model where the baryons are realized as instantons 
would be just as good for our purposes, for example the model of [26] for some 7-brane 
geometries [35], or the AdSs x S 1 model of [32] (the baryons of that model are studied in §6 
of [30]). 

In the baryon-vertex picture, each of the iV c strings connecting the vertex to the flavor 
branes has electric charge 1/Nf under the abelian U(l) subgroup of the U(Nf), so the 
whole baryon has abelian charge N c /Nf. In the instanton picture, the same electric charge 
obtains from the Chern-Simons coupling between the abelian electric field and the non- 
abelian magnetic fields of the instanton. 



ScS ~ w/ ^ ' ^ AT "^ A ) D JfJMi)* g^i > 

n r 

D — - d 5 xA (x) xl(x) (2.26) 

where I(x) is the instanton number density of the magnetic fields. For Nf > 3, the Chern- 
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Simons couplings also endow instantons with non-abelian electric charges. Altogether, 



S cs D fd 5 xA a (x)xQ a J(x) 



(2.27) 



where the net electric charge Q^ has both abelian and non-abelian components; in Nf x Nf 
matrix language, 



for F^ e 
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2 

(2.28) 

Abelian or nonabelian, the like-sign electric charges repel each other; it is this Coulomb 
repulsion between different parts of the same instantons that prevents it from collapsing to 
a point-like D-brane. 

However, the instantons do not grow large because the 5D gauge coupling decreases away 
from the z = hyperplane: a large instanton would spread into regions of space where the 
coupling is weaker, and that would increase the instanton's energy. Instead, the equilibrium 
radius of the instanton scales like 



Oinstanton 



M KK X y/X ' 



(2.29) 



For a holographic model of a baryon, this radius is unrealistically small. Indeed, using the 
p meson's mass as a unit, the real-life baryon radius R^ ~ 3.4M" 1 , while in holography 
a <C Mp 1 ~ M K ^. Moreover, it raises the question of whether we may adequately describe 
so small an instanton using the DBI + CS action, or perhaps we need to include the higher- 
order stringy corrections. 4 On the other hand, assuming the DBI + CS description is OK, the 
radius as in eq. (2.29) allows for consistent expansion of the instanton's parameters in powers 
of A -1 [5]. In particular, the leading contribution to the instanton's mass is 0(X x N c Mk_k) 
while the corrections due to z-dependent gauge coupling and due to Coulomb self-repulsion 
are both 0(N c M KK ) ~ Mj/A. 

To see how that works, consider a static instanton --a time-independent configuration 



4 In string theory, the DBI + CS action for the gauge fields on a D-brane is exact for constant tensions fields 
J- m m however strong. But for the variable tension fields, the DBI action includes all powers of the T mn but 
neglects their derivatives DkJ^mn, DkD v T mni .... It is not clear what effect (if any) such higher-derivative 
terms would have on a small instanton. In a supersymmetric background, the instanton is BPS and its net 
mass is protected against stringy corrections, so the DBI action — or even the Yang-Mills action — gives the 
exact value. But what happens to small instanton in non-supersymmetric backgrounds is an open question. 



22 



of SU(Nf) magnetic fields, plus the electric fields induced by the CS couplings, and the 
scalar fields induced by the tr (<&FF) coupling to the magnetic fields. Since the canonically 
normalized couplings in 5D are 0(A -1 ' 2 ), the leading contribution to the instanton's energy 
comes from the magnetic field. In the DBI approximation 5 

E DB1 = [d 3 x fdz - 2 l x Str (Jdet(K(z)5 fll/ + F fll/ ) - K 2 {z)) (2.30) 

J J *9ym\ z ) V / 

where K(z) = {2 / ko!)~ 1 x gu(z) ~ \j[\M\ K l^). Both #ym(z) and K(z) depend on the z 
coordinate on the distance scale 0(1/Mkk), so for instantons of much smaller size we may 
start with the approximation of constant K and constant 5D gauge coupling. And in this 
approximation, the DBI energy is minimized by the magnetic fields that are exactly self-dual 
(in the 4D space of (x 1 ,^ 2 ,^ 3 , z)); moreover, the DBI energy of an instanton is equal to its 

Yang-Mills energy 

8tt 2 
E-dbi (instanton) = Eym (instanton) = — — (2-31) 

9ym 

regardless of its radius and of the K parameter of the DBI action. 

To determine the equilibrium radius of an instanton we need to calculate its energy to 

the next order in 1/A expansion. To this order, we assume the magnetic fields to be exactly 

self-dual — which allows us to use the YM action instead of DBI even for small instantons 

- but the gauge coupling is z-dependent, and we also account for the electric and the scalar 

fields. For small instantons, we may approximate the 5D gauge coupling as 

= N C \M KK (B + DMl K x z 2 + 0(M^ K z 4 )) (2.32) 



g 2 (z 



for some numerical constants B and D; for the Sakai-Sugimoto model B = C^JTI-ji and 
D = (8£ 3 — 5)/9£ 2 , while other models may have different values. Consequently, the YM 
energy of the nonabelian magnetic fields evaluates to 

£na = N C \M KK x (b + DMl K x Z 2 + DM^ K x y J (2.33) 

where a is the instanton's radius and Z is is the z coordinate of its center. 



In string theory, the complete action for scalar and vector fields living on a D-brane starts with the DBI 
and CS terms but also includes an infinite number of terms involving D n J- mn and the higher derivatives. 
The DBI + CS action is merely a long-distance approximation that allows for strong gauge field tensions T mn 
as long as their gradients are small. It is not clear how well this approximation works for small instantons, 
but at present we do not have a better approximation. 
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The electric potentials Aq couple to the F^F^ products of the magnetic fields while the 
scalar potentials $ a couple to the F llv F ,lv . For the self-dual magnetic fields, both potentials 
couple to the same source /(x, z) x Q^, and the only difference is the coupling strength ratio 
C(z) (cf. eqs. (2.22) and (2.23)). For small instantons we may neglect the z-dependence of 
this ratio and let C(z) ~ C(0) = C, and as long as we do not go very far from the instanton 
(for distances r <C Rkk), we may also neglect the ^-dependence of the gauge coupling. 
Consequently, both electric and scalar potentials become the 4D Coulomb potentials 

$(x, z) = Cx A,(x, z) = *%*«*> /dV [dz> ff^ . (2.34) 

BN C \M KK J J (x - x.') 1 + (z - z'y 

The electric fields lead to repulsive forces between the charges while the scalar forces lead to 
the attractive forces. Altogether, the net Coulomb energy amounts to 



E c = £^N c \M KK x fd 3 xdztr((VA ) 2 - (V$) 2 ) 



(2.35) 



■2W. N c [««,,,, I(X,Z)XI(X.',Z' 



'1 - C z ) x — ^— /d d xdV zddz 



ABXM J (x - x') 2 + (z - z 1 ) 2 ' 

and for a single instanton of radius a it evaluates to 

E c = — — x Nc - . (2.36) 

C B 5\M KK a 2 V ; 

Note: for C < 1 the electric fields are stronger than the scalar fields and the net Coulomb 
energy of the instanton is positive — which makes for a net self-repulsive force that prevents 
the instanton from shrinking to zero radius. In models with C > 1 (assuming they exist), 
the scalar fields would be stronger, the net Coulomb energy would be negative, which means 
a self-attractive force rather than self- repulsive. In such a model, the instanton would shrink 
to zero radius and our approximations would not be valid. 

Altogether, the net energy of an instanton amounts to 

/ fl 2 I _ (J2 

^(instanton) = N C XM KK x [B + DM^ K x Z 2 + £>M| K x 1 x 



2 5B \ 2 Ml K a 2 



(2.37) 

Note that for a 2 ~ 1/(AM^ K ) both radius-dependent terms here are 0(1/ X) corrections to 
the leading term. Minimizing the net energy, we find the equilibrium value of the instanton 
radius 



1 4/2(1 -C 



2^ 



x\7 \" ' , (2.38) 



M KK y/X ' V 5BD 
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the instanton center is in equilibrium at Z = (the bottom of the U-shaped flavor branes), 
and the instanton's mass is 



Mj = N c Mkk x ( \B + \j 2D{l hB ^ + 0(1/A) J . (2.39) 

For the antipodal Sakai-Sugimoto model 

a = </I ^^" 3 ' 2 and M, = N C M KK x ( A + /^ + 0(1/A) ) , (2.40) 
M KK vA \277r V 5 

while other models should have different 0(1) numeric factors. 

We may absorb 2 out of 3 model-dependent parameters B, C, D into a redefinition of the 
A and Mkk parameters of the effective 5D theory, for example 

A ->• -== , M KK -> M = M KK x W - , (2.41) 

VDB S V -D 



thus 



In 2 



9 2 (z) 



N c \M(l + M 2 x z 2 + •••). (2.42) 



For the static instantons or multi-instanton systems we may also get rid of the third model- 
dependent parameter C. Indeed, for the static systems $(x, z) = C x Aq(x.,z), so the only 
effect of the scalar fields is to reduce the net Coulomb force by a constant factor (1 — C 2 ). 
We may simulate this effect without any scalar fields by using different gauge couplings for 
the electric and magnetic fields, g 2 Y — (1 — C 2 ) X g 2 nag , or equivalently by rescaling the time 
dimension x° relative to the space dimensions (rr 1 ,^ 2 ,^ 3 , z), 

t — y , but x — > x, z — > z. (2.43) 

VI -c 2 
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Consequently, the static instanton's energy becomes 



E nct = N C \M + N C XM 3 x d 3 xdzl(x,z) x z 



+ TTT7 x /^ 3 x dz rfV dz . ^' *(* ; ,Z { + subleading (2.44) 

4AM J (x - x') 2 + (z - z') 2 v ; 

N C \M + N C XM 3 x ( Z 2 + — ) H — c — + subleading, 

V 2 / 5AMcr 



hence in equilibrium Z = and 



^2/5 
Mv 7 ! ' 



(2.45) 



Besides the radius a and Z, the instanton has other moduli — the X 1 ' 2,3 coordinates of 
the center (which corresponds to the baryon's coordinates in 3D) and 4Nf — 5 orientation 
moduli in the SU(Nf) gauge algebra; and the net energy is degenerate in these moduli to all 
orders in 1/A. For finite N c — even if it's very large — one should quantize the motion of the 
instanton in those moduli directions. Consequently, a holographic baryon acquires definite 
spin and isospin quantum numbers; for Nf = 2 the baryons have / = J = 0, 1, 2, . . . , 4^ for 
even N c or I = J = |, |, |, . . . , 4p for odd N c [5, 30], and there are similar spin -h- flavor 
multiplet relations for Nf > 2. 

However, in multi-baryon systems interactions between the baryons break the rotational 
and flavor symmetries of individual baryons, and only the overall 5*0(3) and SU(Nf) sym- 
metries remain unbroken. In holography, the multi-instanton systems suffer from the same 
problem: the magnetic fields of multiple instantons interfere with each other, which spoils 
the degeneracy of the net energy with respect to orientations moduli of the individual instan- 
tons. In the large N c limit, this effect becomes more important than the quantum motion in 
the moduli space. 

Consequently, in our research program into holographic nuclear matter and related issues, 
we stick to classical static instantons with definite classical orientations in space and in 
SU(Nf). From the quantum point of view, such instantons are superpositions of states with 
different spins and isospins (or rather 5*0(4) and SU(Nj) quantum numbers). We do not 
even care is the instantons are bosons or fermions — at our level of analysis it simply does not 
matter. Minimizing the classical energy of a classical multi-instanton system with respect 
to classical positions and orientations of all the instantons is already a very hard problem - 
see the next section for the highlights, — and we do not wish to delve into the quantization 
of moduli until we have the classical issues under control. 
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2.5 Multi— Baryon Systems 

In the large N c limit, nuclear forces between the baryons are dominated by the static po- 
tentials. Holographically, a static system of A baryons corresponds to a time-independent 
configuration of the non-abelian magnetic flavor fields F" (x, z) (n,v = 1,2,3, z) of net 
instanton number A, 



d xdz 
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tr(F KA F ; 



k\" ixv) 



A = ^baryons, 



(2.46) 



accompanied by the Coulomb electric Aq(xl,z) and scalar $ a (x, z) potentials induced by 
their Chern-Simons and QFF couplings to the magnetic fields [2]. The whole configuration 
should minimize the net DBI + CS energy of the system subject to the constraint (2.46). 

In the A — > oo limit, the DBI energy of the magnetic fields is 0(A) while the net effect 
of the electric and scalar fields is only 0(1). Moreover, the magnetic fields are concentrated 
within 0(A -1 / 2 ) distance from the X 4 = hyperplane, so to the leading order we may 
approximate the 4 + 1 dimensional spacetime as flat. Thus, 
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(2.47) 



where we neglect the z dependence of the g Y M anc ^ K = ^n (^)/(27ra; / ). Similar to the 
Yang-Mills energy of an A instanton system, this leading-order DBI energy is minimized 
by the self-dual configurations of the magnetic fields F" (x, z). In fact, all such self-dual 
configurations (of the same instanton number A) have the same leading-order energies 
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(2.48) 



The self-dual configurations form a continuous family parameterized by A x ANf moduli, 
which correspond to the locations, radii, and SU(Nf) orientations of the A instantons. But 
the leading-order energy (2.48) does not depend on any of these moduli. 

Fortunately, the sub-leading corrections to the net energy lift the degeneracy of the 
leading order, which provides for the 0(A°) interactions between the baryons. To work 
out such interactions we need the degenerate perturbation theory for the magnetic field 
configurations and their energies. At first order of this perturbation theory, we (a) limit the 
F" (x, z) configurations to the degenerate minima of the leading-order energy function, i.e., 
to the self-dual magnetic fields; (b) calculate the 0(A°) corrections AE for the energies of 



27 



these configurations; (c) minimize the AE among the self-dual configurations. At the next 
order, we would calculate the 0(\~ l ) corrections AF" (x, z) to the magnetic fields — which 
would no longer be self-dual -- and then use such corrections to calculate the net energy 
to the order 0(A _1 ). But in this article, we limit our calculations to the 0(X ) interaction 
energies and zeroth-order magnetic fields, so we shall stop after step (c). 

In more detail, to get 0(A°) interactions between A baryons, we proceed as follows: 

1. General self-dual magnetic field configurations obtain via ADHM construction [11] in 
terms of A x A and A x Nf matrices obeying certain quadratic constraints. Our first 
task is to solve these constraints and write down the ADHM matrices in terms of the 
instantons' locations, radii, and orientations. 

2. Given the ADHM matrices, we work out the instanton number density profile 



16^ 



I(x,z) = -^tr(F KX F^), (2.49) 



and for Nf > 2 also the non-abelian adjoint density 

r^x) = j^xd abc F b KX F^. (2.50) 

3. Next, we calculate the 0(A° x N C M) corrections to the net energy of the system. There 
are two sources of such corrections: the z-dependence of the 5D gauge coupling, and 
the Coulomb electric and scalar potentials induced by the Chern-Simons and the &FF 
couplings to the magnetic fields, thus 

AE = AE NA + AE C . (2.51) 

The ^-dependent 5D gauge coupling changes the DBI energy of the magnetic fields by 

A£ NA = N C M x d 3 xdz\M 2 z 2 xl(x,z), (2.52) 



while the Coulomb energy depends on the Nf. For Nf = 2, the U(2) Chern-Simons 
and &FF terms couple the SU(2) magnetic fields to the U(l) Coulomb fields only. 
Consequently, the Aq and the $ a fields are abelian and couple to the instanton density 
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J(x, z), and their net energy is simply 4 + 1 dimensional Coulomb energy 

N C M /"„ , /•- , /(xi,^i) x i"(x 2 ,z 2 ) 

A£ c = -^— x /rf^rfx! d 3 x 2 dz 2 \}' V ±j^ ^-. 2.53 

A J J (AM 2 ) x ((xi - x 2 ) 2 + ( Zl - z 2 f) 

For JV/ > 2, the £/(iV F ) C/(2) Chern-Simons and $FF terms couple the SU(N F ) 
magnetic fields to both abelian and non-abelian electric and scalar fields; the abelian 
Coulomb fields are sourced by the instanton density J(x, z) while the non-abelian fields 
are sourced by the adjoint density I a (x, z). Altogether, the net energy of these Coulomb 
fields is 

N r M f „ f , jfcJ(xi,2i) x 7(x 2 ,^ 2 ) + 4/°(xi,«i) x J a (x 2 ,z 2 ) 



A£ c = — ^ — x d 3 x 1 dz 1 d 3 x 2 dz 2 — 



(AM2)x((x 1 -x 2 ) 2 + (^ 1 -^ 2 )2) 

(2.54) 



4. Steps 1, 2, and 3, give us AE as a function of baryons' locations, radii, and SU(Nf) 
orientations. Now we need to minimize the AE with respect to all these moduli. 

This 4-step procedure is fairly straightforward for a few baryons - - cf. calculations of 
the 2-body nuclear forces by Kim, Sin and Zahed [33] and by Hashimoto et al [34] , - - but 
it becomes prohibitively difficult for large numbers of baryons and outright impossible for 
infinite baryon crystals. At best, we can survey a small subspace of the A- instanton moduli 
space and try to minimize the AE over that subspace. For example, we may assume that 
all the instantons have the same radius, that their centers form a periodic lattice of some 
particular symmetry, and that the orientations of the instantons also form some kind of a 
periodic pattern; this gives us an ansatz for all the A x 4Nf moduli in terms of just a few 
overall parameters, and we can try to calculate and minimize the AE as a function of these 
parameters. However, any such ansatz is likely to miss the true lowest-energy configuration 
of the system. Indeed, in condensed matter guessing the crystalline symmetry of some 
substance from the properties of the individual atoms is a game of chance with poor odds, 
and there is no reason why the instanton crystals should be any simpler. Moreover, even if we 
could somehow guess all the symmetries of the instanton crystal, actually working through 
steps 1-4 is impossible without additional approximations (besides A>1). For example, in 
our previous paper [1] we needed the (lattice spacing) ^> (instanton radius) approximation 
just to solve the ADHM constraints for the zigzag-shaped chain of instantons, and even for 
a straight chain we needed the spacing ^> radius approximation to calculate the Coulomb 
energy AEq- 
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In this article we explore a shortcut around steps 1, 2, and 3. In section 3 we shall see that 
when the distances between the instantons are much larger than their radii, the interactions 
between the baryons are dominated by the two-body forces, 

interactions » g S F * ( X ™ ~ X ™> orientation m, orientation n ) (2.55) 

m,n=l,...,A 

for a manageably simple function F 2 of the two instanton's positions and orientations. Con- 
sequently, we can minimize the net interaction energy over the entire moduli space of the 
multi-instanton system using a simple numerical simulation: Starting from a random set of 
instanton positions and orientations, use the steepest descent algorithm to find the nearest 
local minimum of the net energy; repeat this procedure for different random starting points 
to find other local minima; eventually, we find all the local minima, compare their energies, 
and identify the global minimum. In our next paper [41] we shall use this algorithm to 
construct the 2D and 3D instanton lattices starting from scratch, i. e., from the F 2 function. 
In the present article, we focus on the ID lattice (§4) and on the first step in the transition 
form ID to 2D -- the zigzag (§5). In both cases, we assume the lattice geometry for the 
instanton positions, but we use the numeric simulation to find the lowest-energy pattern of 
their orientations. Once we know all the patterns that show up, we use them as ansatz's 
(with a few parameters) for which we calculate the net energy as analytic functions of the 
parameters. In this way, we map the boundaries between different orientation patterns much 
more accurately than we could do it in the numerical simulation. 

To make the baryons form a ID lattice instead of spreading out in three dimensions, 
we curve the x 2 and X3 dimensions of the flavor brane similar to the curvature of the z 
coordinate. In terms of the effective 4 + 1 dimensional theory, this corresponds to the 5D 
flavor gauge coupling depending on the x 2 and the X3 as well as the £4 = z, 

-^— = N C \M x (l + M\x\ + Mlx\ + M\x\ + 0(M 4 x 4 )). (1.2) 

#5 K x ) ^ ' 

This gauge coupling acts as a harmonic potential for the instantons which pulls them towards 
the x\ axis, so at low densities the instantons form a ID lattice along the x±. At higher 
densities, the instantons push each other away from the x\ axis and form more complicated 
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2D or 3D lattices, starting with the zigzag 

„2 



*x A' 



■/*<• ■•>•••■ 

x^; x^r. ... v y " xl 



To make sure the transition from the straight chain to the zigzag happens for lattice spacings 
much larger than the instanton radius (which is required by the two-body force approxima- 
tion), we assume M 2 <C M 3 , M4. 

This setting is similar to what we have used in our previous paper [1], except for a 
change of notations 6 and allowing for M3 7^ M4. In §5 we shall see that M3 7^ M4 makes a 
big difference for the instanton orientation patterns in a zigzag: the lowest-energy patterns 
are abelian for M3 m M4 but non-abelian for anisotropic M3 < O.9M4. Even the order 
of the phase transition from the straight chain to the zigzag depends on the M3/M4 ratio: 
second-order for ratios > 0.63 but first-order for ratios < 0.63. 



3 Two-Baryon versus Multi-Baryon Interactions 

In real-life nuclear physics, besides the two-body nuclear forces due to meson exchanges, 
there are significant three-body forces, and presumably also four-body forces, etc., 

A 

Nucleus = ]T£ lbody H + i J2 H 2hod y(m,n) + i J2 H 3hod y(£,m,n) + ••■ 

n=l different different 

m,n=l,...A £,m,n=l,...A 

(3.1) 
where n stands for the quantum numbers of the n th nucleon. Likewise, in the holographic 
nuclear physics interactions between multiple baryons include two-body forces and also three- 
body, four-body, etc., forces. Even in the classical infinite-mass limit where the holographic 



'Back in [1], we had 

-^— = 7V C AM x (l + M' 2 (xj +x 2 2 ) + M 2 xl + 0(MV)Y M' > M, 



M ncw 



so the zigzag was in the {x%, a; 4 ) plane, with the long axis along the X4. Translating between the old notations 
and those of this paper, we have ( Xl ,x 2 ,x 3l x 4 ) old = (x 4 , x 3 ,x 2 ,x 1 ) ncw , M old = M 2 ncw , M' old = Af 3 ncw = 

nc 
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baryons become static instantons of the SU(Nf) gauge fields in 4+1 dimensions, the potential 
energy (due to g^ 7^ const and due to Chern-Simons interactions) of an A-instanton system 
has form 

S net (l,2,...,A) 
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with significant three-body, etc., terms. 

What about the relative magnitudes of the 2-body, 3-body, 4-body, etc., interaction 
terms? When the baryons are packed cheek-by-jowl till their instanton cores overlap and 
merge, we expect all the n-body forces to have comparable strengths. But in the opposite 
low- density regime of baryons separated by distances much larger then their radii, the two- 
body forces dominate the interactions, while the multi-body forces are smaller by powers of 
(radius / distance) 2 . The purpose of this section is to prove this fact. 

Before we delve into the technical issues, let us briefly summarize our argument. First, 
consider the simplest case of Nf = 2 and intermediate-range distances between the baryons, 

< \X m -X n \ < — , (3.3) 



MVxJ ' m n ' M 

which allows us to treat the 5D holographic space as approximately flat. For Nf = 2 
the holographic baryons are instantons of the £77(2) magnetic fields, which source the U(l) 
electric and scalar fields via CS and tr^F^F^) couplings. When those instantons are small 
and separated from each other by large distances \X m — X n \ ~ D ^> a, their interactions 
come from two sources: 

1. Direct Coulomb repulsion (electric — scalar) between nearly-point-like abelian charges in 
4 + 1 dimensions, 

^direct N c V^ 1 /q A) 

which is a manifestly two-body interaction. 
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2. The interference between the instantons changes the distribution of the instanton number 
density in space, 



32^ 

A 



T ( x ) = i^tr(F^(x)F pa (x)), (3.5) 



\A instantons / \ 



J^ /standalone^) + ^interference (^ ^ 



which in turn changes the self-interaction energy (Coulomb and non-abelian) of each instan- 
ton by an amount comparable to (3.4). Within an instanton — i. e., at 0(a) distance from 
some instanton's center X% — the interference from the other instantons should be relatively 
weak, 



a 
D 



2 



A ^interference / \ w rstandalonc / \ /o t\ 

AI (x) ~ — x I n (x), (3.7) 



so there should be some kind of a perturbation theory for it. At the first order of such per- 
turbation theory, the AI arises from interference between the un-perturbed standalone-like 
instantons, so we expect it to be a sum of pair-wise interferences from the other instantons, 

A/ lstordcr (xnearX n ) = /f dalmc (x)x ^ jj(i) (3.8) 



a 2 



H%(x) ~ O ( - ^- 2 ) (3.9) 



where 

F« (t.) ~ (1 ( 

\v Y 

and depends only on the instantons j^n and j^m - - i.e., on their positions, radii, and 
orientations -- but not on any other instantons. At the second-order, we expect to include 
the interference between the first-order AI and the additional instantons, so at this order 
we obtain 3-body effects, 

A/ 2ndordcr (:rnearX„) = if andatoM ( x ) x ^ j£}, m (x) (3-10) 

but 

^ (2) ~ £j « ^. (3.11) 

Likewise, the higher orders may involve more and more instantons, but the magnitudes of 
such high-order interference effects are suppressed by the higher powers of (a 2 /D 2 ). 
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Now consider the Coulomb self-interaction of the instanton #n, 

N r 1 



^" ( "> = 4AS X 7 n < 3 ' 12 ' 



where p n is the instanton's effective charge radius, 



\ = ff fr^pM&L. (3.13) 

Pi J J \xi-x 2 \ 2 

xi, X2 near V n 



A standalone instanton has 

1 4/5 



Pn °n 



(3.14) 



but interference from the other instantons should change this radius by small amount of 
similar relative magnitude to A///;f andalone , thus 

14/5**1 / 

-if -> -jr + A «> A « ~ H- ( 3 - 15 ) 

Pn u n -^ 

which changes the instanton's Coulomb self-interaction energy by 

A£ s r eli (n) = — xA n ^— . (3.16) 

c v ; 4AM " AMD 2 v ; 

Note that this effect has a similar magnitude to the direct Coulomb repulsion (3.4) between 
the instances. 



Moreover, the charge radius correction A n is linear in the A7 m terterence at x near the X, 



n ■ 



hence in light of eq. (3.8), the leading-order contribution to the A n is a sum of pair-wise 
interferences from the other instantons, thus 

An = J2 A nl + 0(a 2 AD 4 ) (3.17) 

where each A n , m depends only on the instantons #n and #ra. Consequently, the leading 
effect of the interference on the net Coulomb self-energy of all the instantons has form 

A interference£seif = £^ clf [interfering] - ^[standalone] 



4AM 

n,m=l,...,A 
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V A« + O (— x —} (3-18) 

, — 1 A \ / 



where the leading terms act as two-body interactions between the instantons. 

Besides the Coulomb energy, the non-abelian energy is also affected by the interference 
between the instantons, 



interference 



E NA = N c \M 3 x fd 4 x(x 4 ) 2 xAI mtcllcicncc (x), (3.19) 

but this time the integral should be taken over the whole 4D space, including both the 
instantons and the inter-instanton space. Indeed, we shall see later in this section that 

at 0(a) distance from an X n , AI 



a f (3.20) 

at 0(D) distances from all the X m , AI ~ — — , 

and both kinds of places make 0(a 4 /D 2 ) contributions to the integral (3.19). Moreover, in 
both kinds of places, the leading terms in the AI(x) is a sum of independent two-instanton 
interference terms, 

^interference^) = 1^2gi(x) + Subleadhlg. (3.21) 

Although our heuristic argument for such decomposition near instanton centers does not work 
in the inter-instanton space, we shall prove later in this section that the decomposition (3.21) 
works everywhere in the 4D space. Therefore, the non-abelian interactions between the 
instantons due to interference are also dominated by the two-body terms 

interference^ = 1 ^ JV C AM 3 fd 4 X (x 4 ) 2 X I n>m (x) + Subleadhlg. (3.22) 

Thus far we have assumed Nf = 2. For larger numbers of flavors, interference between 
magnetic fields of the SU(Nf) instantons is harder to describe mathematically, but its qual- 
itative features for a ^ D remain exactly as in eqs. (3.5) through (3.21). In particular, 
the effect of interference between small instantons separated by large distances on the den- 
sity profile I(x) is dominated by the two-instanton interference terms as in eqs. (3.21) and 
(3.8-11). Consequently, the instantons' non-abelian and Coulomb interaction energies suffer 
small corrections from the interference that give rise to predominantly-two-body interactions 
between the instantons as in eqs. (3.18) and (3.22). 



35 



The qualitatively new aspect of Nf > 3 comes not from the SU(Nf) self-dual magnetic 
fields of the multiple instantons but from the electric and scalar Coulomb fields induced by 
the magnetic fields via Chern-Simons and tr(^F' muu F^ v ) interactions. For Nf = 2 such 
electric and scalar fields are purely abelian, A , $ e U(l) C U(2), but for Nf > 3 they have 
both £7(1) and SU(Nf) components. For a single instanton, all fields — magnetic, electric, 
and scalar - - belong to the same Cartesian U{2) subgroup of the U(Nf) - - in Nf x Nf 
matrix language, 
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(3.23) 
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- so the Coulomb self-interaction of a standalone instanton works exactly as for Nf = 2. 
On the other hand, Coulomb interactions between multiple instantons whose magnetic fields 
span different 577(2) subgroups of the SU(Nf) become much more complicated. 

Moreover, in a general multi-instanton background, the Coulomb fields A°(x) and $(x) 
do not behave as simply 1/r 2 . Instead, they satisfy gauge-covariant equations for the matrix- 
valued fields, 
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(3.24) 



and ditto for the scalar $(x). Consequently, the Coulomb force between two sources (such 
as instantons) depends on the magnetic fields in the space between them. In other words, for 
Nf > 3 the direct Coulomb interactions are no longer two-body; instead, the force between 
two instantons depends on the other instantons in the system. 

Fortunately, for small instantons separated by large distances, the magnetic fields are 
concentrated in small volumes of 4-space, so their effect on the Coulomb field propagation 
becomes perturbatively weak: 
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(3.25) 



+ 0(a J D ) two-instanton terms + 



where b and c are adjoint indices of the U(Nf). The leading term here is the ordinary 
Coulomb propagator in 4D space, thus in the a <^ D limit, the direct Coulomb interactions 
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between small instantons are predominantly two-body, 

^direct _ N c y^ overlap (n,m) 

b ~ 4AM 4- |X£-X&|* [6 - Ib) 

where overlap(n, m) > is the overlap between the S77(2)'s spanned by the magnetic fields 
of the instantons j^-n and j^m. 

The bottom line is, for all Nf > 2, all the interactions between small instantons separated 
by large distances are predominantly two-body. However, for Nf > 3 the dependence of such 
two-body interactions on the instantons' orientations becomes much more complicated than 
for N f = 2. 

Finally, let us relax our second un-necessary assumption of intermediate-range distances 
between the baryons and consider what happens at longer distances D ~ 1/M. In this 
regime, the curvature of the fourth space dimension - - and especially the x 4 -dependence 
of the 5D gauge coupling - - can no longer be treated as a perturbation. Consequently, 
the magnetic fields of an instanton or a multi-instanton system are no longer self dual in 
the inter-instanton space - - although inside the instanton cores or at intermediate-range 
distances r <C 1/M from instantons they remain approximately self-dual. Therefore, the 
interference between very distant instantons is no longer governed by the self-dual ADHM 
solutions. Instead, we must work it out the hard way: figure out how the magnetic fields of 
an instanton propagate through the curved 4 + 1 dimensions towards the other instantons, 
and then find out how such fields disturb the other instantons' cores. 

Fortunately, we do not need a hard calculation to see that at large distances from an 
instanton its magnetic fields are very weak. Indeed, even in flat 5D space the fields weaken 
with distance as A a ~ a 2 /r 3 (in the IR-safe singular gauge), so at r > a they are so weak 
that the field equations become effectively linear. In the curved space, we may decompose 
the weak 5D gauge fields into 4D mesonic fields, hence at large distances r > 1/M from an 
instanton, its fields become Yukawa-like 

2 2 

A«(r,Z) ~ ^xJ2^k(z)^ k (z- mst )e- m ^ r « ^ . (3.27) 

k 

Consequently, at the location X% of any particular instanton, the background fields from the 
other instantons are very weak, and their effect on the instanton #n itself can be adequately 
accounted by the first-order perturbation theory. In other words, the effects of other instan- 
tons j^m 7^ n on the instanton j^n are weak and add up linearly! For the self-interaction 
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energy of the n th instanton, this means 

£{n) « ^(standalone) + ^ A m S(n) (3.28) 

where the second term gives rise to two-body interaction energies 

interference («, "0 = ^(n) + A n £(m). (3.29) 

We expect the two-body terms to be rather small — in addition to the usual 5D l/\X m — X n \ 2 
factors they should carry Yukawa-like exponentials e~ mr (or rather sums of such exponen- 
tials), but these are the leading interactions due to interference. To three-body or multi- 
body interactions follow from higher-order perturbation by very weak fields (3.27), so they 
are much smaller that the two-body interactions. 

As to the direct Coulomb interactions between the holographic baryons via electric or 
scalar fields, at large distances \X m — X n \ ~ 1/M they also decompose into sums of Yukawa 
forces. Moreover, since the scalar mesons generally have different masses from the vector 
mesons, the attractive potential due to scalars may have a different r dependence from the 
repulsive potential due to vectors. Thus, for a right model, the net two-body force between 
two holographic baryons may become attractive at large enough distances between them. 
But regardless of the model, the direct Coulomb interactions are always manifestly two-body 
for Nf = 2, while for Nf > 3 the multi-body terms exist but become very small at large 
distances between the baryons. Thus, 

£direet = \ J2 °verlap(m, n) x V {2 \\X m - X n \) (3.30) 

where the precise form of the potential V^ 2 \r) is model-dependent, but the two-body form 
of the direct interactions is quite universal. 

* * * 

At this point we have finished our heuristic arguments and about to start a more technical 
analysis of the interactions between the holographic baryons. To keep our presentation 
relatively simple, we assume just two quark flavors and intermediate-range distances between 
the holographic baryons — D ^> a but D <C 1/M . 

For Nf = 2 the non-abelian gauge symmetry is £77(2), so the ADHM data for an A- 
instanton configuration consists of a quaternionic array Y n (n = 1, . . . , A) and a symmetric 
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quaternionic A x A matrix T mn . Equivalently, we may use A real instanton radii a n = 
\Y n \, A SU(2) matrices y n (equivalent to unimodular quaternions Y n /a n ) parameterizing 
the instanton's orientations, and 4 real symmetric matrices r^ n whose diagonal elements 
T^ n = X% are the locations of the instanton's centers. Finally, the off-diagonal matrix 
elements T^ , = <x^ nn follow from the other parameters by solving the ADHM equations 

S(( rtr L„ + ^n) = ( 3 - 31 ) 

or equivalently 

<4 rM > r 1 m n + °- a « x trfe„(-^ a )) = (3.32) 

where 77" is the 't Hooft's tensor mapping between the 5 , f/(2) gauge and the SU(2) L inside 
Spin(4) = SU(2) L x SU(2) R , 

*& = 0, V« = -%, ^ = +^ a , vtj = e aiJ , a,i,j = 1,2,3. (3.33) 

The ADHM data are somewhat redundant --an 0(A) symmetry acting on all the T^ n 
and Y n = a n y n does not change any physical properties of the multi-instanton data. This 
symmetry includes Z 2 which flips the instanton orientations y n — ¥ —y n (independently for 
each n). It also includes small SO (A) rotations that change the off-diagonal elements by 
boi mn = e mn (X m — X n y + 0(e 2 ). To eliminate these rotations, the ADHM equations (3.32) 
for the off-diagonal elements should be combined with additional constrains (one for each 
m 7^ n), for example 

Vm^n: (X m - X n )^< n = 0. (3.34) 

For large distances \X m — X n \ ~ D between the instantons, D ^> a n , we may solve the 
ADHM equations (and the constraints (3.34) as a power series in a 2 /D 2 : 

u _ -pM _ (1) , (2) 1 (3) 1 . . . (Q QZ.^) 



a 



gL = ^^ ^ x \a m a n tr(yly n {-ir a )) = 0(a 2 /D), (3.35b) 



iimn I \r y 

Am — A; 



(2) VuvV^-rn ^n)v sr^ a (1) (1) , 4 , 3, , . 

%L = — , y _ Y |2 x ^ ri K x®J m a xL = °( a ID 6 ), (3.35c) 

z?=m,n 

«S» = - ^\f m ~vt >< £ ^«2L«S. = ^(a 6 /^ 5 ), (3.35d) 

-^ 7TJ. -^ TI 



£^m,n 
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Note that for each off-diagonal matrix element, the leading term ajxmn in this expansion 
depends only on the instantons #m and #n (i. e. , on the positions, radii, and orientations 
of only these two instantons) while the subleading terms a^mn, aji™, • • • involve additional 
instantons. 

Given the ADHM data T^ n , a n , and y n of a multi-instanton system, the 4D instanton 
number density (3.5) can be obtained as 

I(x) = -— DDlogdet(L(:r)) (3.36) 

where L(x) is a real A x A symmetric matrix 

L mn {x) = ^(r^ m -x At ^ m )(r^-x M ^ n ) + ja m a n tr(yl n y n ). (3.37) 

e 

Thanks to the double D'Alambertian in eq. (3.36), several moments of the instanton density 
may be obtained via integrating by parts: 

d 4 xl{x) = A, (3.38a) 

fd 4 xl(x)xx u = tr(T u ), (3.38b) 

fd 4 xl(x)xx^x u = tr(FT^) + |^tr(T) (3.38c) 

where T mn = \a m a n tr(yj„y n ), (3.38d) 

ki(i)xAv = |(r A {r^,r}) + ^tr(rr) (3.38e) 

+ ±5 Xu tr (T»T) + |^tr(r A T). 

The quadratic moment (3.38c) for \x — v — 4 is particularly important since it gives us an 
exact formula for the non-abelian energy of the system directly in terms of the ADHM data, 

E NA = N C \M 3 x fd 4 xl(x) x (x 4 ) 2 = N C \M 3 x (tr(r 4 r 4 ) + ±tr(T)) 

= NAM' X (j2{( T nn) 2 + Fun) + Y.Kn)"\ (3-39) 

\i=n rriy^n / 

= N c \M 3 J2((Xnf + K) + N C XM 3 J2Krf- 



m^n 
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Obviously, the first sum on the last line here is the sum of individual instantons' potential 
energies due to their radii and locations (relative to the x A = hyperplane) while the second 
sum comprises the interactions between the instantons, 

^interaction = N c \ M 3 ^ (a^f . (3-40) 

This formula is general and applies to any multi-instanton configuration, the problem is cal- 
culating the off-diagonal matrix elements a^ n = r^ , . Fortunately, for the small instantons 
separated by large distances 0(D) ^> a, those off-diagonal matrix elements are given by the 
expansion (3.35.a-d) in powers of (a 2 / ' D 2 ). Moreover, for each m ^ n the leading term cv^mn 
in this expansion depends only on the positions and orientations of the instantons j^m and 
j^-n but does not depend on any other instantons. Thus, to the leading order in a 2 / D 2 , the 
non-abelian interactions (3.40) between the instantons are two-body, i.e. involve only two 
instantons at a time: 

fraction w 1_ £ S ^(m, n) ~ iV c AM 3 -^ (3.41) 

while the multi-body interactions are 0(N c XM 3 a 6 / D A ) or weaker. Specifically, 



S 2 N h f y (m,n) = 2N C \M 3 x ( a m a n x ^ Xm _ * w) " x \ tr(yly n (-ir a )) 



2 



NAM 3 a 2 a 



3 - 2 a 2 



0\Y _ "V^ 1 2 

Zj \ yv m A n|4() 



x tr 2 (yly n (-iN mn • r) 3D ) (3.42) 



where N mn is the 3- vector part of the unit 4- vector 



N^ n ee {N mn ,N^ n ) = g _ y . (3.43) 



Now consider the Coulomb energy 



4AM J J \xi-x 2 \ 



£c = iih lld^d Wr^ . (3.44) 



of the multi-instanton system. Alas, we cannot obtain this Coulomb energy directly from 
the ADHM data via integration by parts or some other clever trick but have to calculate 
it the hard way: Derive the instanton density profile from eq. (3.36), decompose I(x) into 
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individual instantons plus interference terms, and finally perform the integral (3.44) for some 
approximation to the integrand. 

Our starting point is the L mn (x) matrix (3.37) which governs the instanton density profile. 
For a system of small instantons separated by large distances, the diagonal elements of this 
matrix are much larger than the off-diagonal elements. Indeed, the diagonal matrix elements 
are given by 

L nn (x) = ({T»-x»f) nn + T nn = {X»-x»f +J2K1 2 + °n ( 3 - 45 ) 

and their magnitudes are generally 0(D 2 ). To be precise, in the inter- instanton space all 
of the L nn (x) are 0(D 2 ) or larger, while in the neighborhood of an instanton #m the 
corresponding L mm (x) may become as small as 0(a 2 ) but all the other L nn (x) remain 0(D 2 ). 
On the other hand, the off-diagonal matrix elements are only 0(a 2 ) <C D 2 , 

J-'m^in — \\*- % ) ) mn ' -'-mn 

= (X m + X n - 2x)» < n + J2 a im a in + a m a n x\ti{yly n ) (3.46) 

= 0(a 2 ). 

This hierarchy between the diagonal and the off-diagonal matrix elements of L allows us to 
evaluate logdet(L) as a power series in in the ratios of the off-diagonal to diagonal elements: 



(3.47) 



logdet(L) = log det (L diag ) + logdet(l + L diag L ffdia g ) 
= trlog(L diag ) + trlog(l + L^l g L oSdiSig ) 
= trlog(L diag ) + tr(L d 4 g L offdiag ) - | tr(L d 4 g L ofMiag L d 4 g L ofMiag ) 

+ 3 tr (-^di a g-^offdiag-f' d i ag -^offdiag-^ d i ag -f'offdiag) + ••• 

Ei / t \ r> 1 \ ^ * J mn* J nm l \ r *- 'fan*- 'mn-'-' 'n£ 
log(L nn ) + - -^ + _ ^ + . . . . 

/ - u mm- u nn ,.„ , - u U- u mm- u nn 

n m^tn dirtcrcnt 

£,m,n 

Consequently, the net instanton density profile of the system may written as a sum of one- 
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instanton, two-instanton, three-instanton, etc., terms: 

-1 



I x) 



DD logdet(L) 



16tt 2 

/ j -'-n \ X ) ~r 2 / j ran 



X 



+ i V x 1 

+ 6 Z^ X ^ 



tmn\ X ) 



m^n 



different 
Lm.n 



(3.48) 



where 



X«, 



1 



x 



X (2) (x) 



■^lran\ X ) 

etc., etc. 



167T 
+ 1 

16vr 2 

-2 

16tt 2 



2 DDlog(L nn ( a ;)) ; 



DD 



DD 



l-'mn-'-'nm 

1 T ' 

- LJ mm- L -'nn 
J-'£m J-'mn J-'nt 

l-'U -L-'mm l-'nn 



(3.49) 
(3.50) 
(3.51) 



The one-instanton terms (3.49) here have the form of stand-alone instanton profiles. Indeed, 
rewriting eq. (3.45) as 



x 



2 ~2 

T — Q, 



where r 2 = \x^ — Xi"\ 2 and al 



Via" I 2 

rriy^n 



(3.52) 



a 2 + 0(a 4 /X> 2 ), 



we obtain 



-1 
16^ 



DD log(L r 



16tt 2 



d 2 3d 

+ 



()r- 



r„ dr. 



log(r 2 + a 2 ) 



(6/tt 2 ) a 4 
(r 2 + a\ f 



(3.53) 



which is precisely the density profile of a standalone instanton of radius a n . Note that this 
radius is slightly larger than the original instanton radius a n — this is a hidden effect of the 
interference from the other instantons m ^ n — while the visible effects of the interference 
come via the two-instanton terms Zmn(x), the three-instanton terms X^J n (x), etc. Deriving 
explicit formulae for all these terms is a straightforward (albeit rather tedious) exercise, 
and the results are too cumbersome to present here. Instead, let us simply estimate the 
magnitudes of the multi-instanton terms. 

In the inter-instanton space — i. e. at 0(D) distances from all the instanton centers - 
we have 



^-'mnJ-'nm 
J-'mm^-'nn 



LjPrn i-Jm.n i-^7 



D 4 ' L Pf L 



~ ^t , etc. (3.54) 



U-'-'mm-'-'nn 



£>6' 
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Moreover, all these expressions change rather slowly with x — on the scale x ~ D — so in 
the context of such expressions d^ = 0(1/ D). Consequently, 

□ i — i J- J rnnJ-'nm 0, — — ±j£ m ±j mn ±j n £ d /„__-, 

D 7 — 1 75s' UU T1 — t 75^' •"' (3 ' 55) 

- u mm lj nn LJ ±J U 1J mm- u nn J - / 

and therefore 

4 6 

In the inter-instanton space 1$ n {x) ~ -=-g , lfm n {x) ~ j^ , (3.56) 

As promised in eq. (3.21), the interference between small well-separated instantons is dom- 
inated by the two-instanton terms. In fact, in the inter-instanton space the two-instanton 
terms have similar magnitudes to the one-instanton terms, 

ZffiO*) ~ ^ ~ ^Hx). (3-57) 

However, the entire I(x) in the inter-instanton space makes a negligible contribution to the 
net Coulomb energy of the system, 

S c [llS] = fd 4 xl(x) x -^V- l —— ~ D 4 x —^ x — ^- ~ -^— R - (3.58) 

CL J 7 KJ 2\M ^ \x» - X%\ 2 D s XMD 2 AMD 6 V ; 

us n 

Therefore, the net Coulomb energy of the multi-instanton system may be summarized as 



£, 



, AT, 

c 



net _ J ' c 



E i + E |v^ Y M|2 + °(« 4 /^ 6 )) , (3-59) 



4AM 1 ^ p? ^ |X& - X s 



where the first term inside the '{•••}' comprises the self- interactions of individual instantons 
- but with the effective charge radii p n that include the effects of interference from the 
other instantons, — the second term is the direct Coulomb repulsion (3.4) between different 
instantons approximated as point charges (or equivalently, as compact spherically-symmetric 
charges), and the third term accounts for the 0(a 4 /D 8 ) instanton density in the inter- 
instanton space. Note that the self-interaction terms dominate the net Coulomb energy, so 
even relatively small corrections to the effective charge radii pi can have an effect comparable 
to the direct Coulomb repulsion between the instantons. 

To work out the effects of interference on the charge radius p n of the n th instanton we need 
to estimate (and then evaluate) the interference terms X)J, X kim) etc., for x^ near the n th 
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instanton's center Xfi. In this neighborhood, the estimates (3.56) do not work — or rather 
they work only for the interference terms between the other instantons #k, £, m, . . . 7^ n. 

But the terms iLl, 2^J n , etc., that do involve the n th instanton itself turn out to be 

much larger than (3.56) because of the 1/L nn (x) factor that happens to be 0(l/a 2 ) instead 

of 0(1/ D 2 ). Moreover, this factor depends on x much more rapidly than all the other factors 

- on the scale of x ~ a instead of x ~ D -- so taking the space derivatives of this factor 

produces much larger results than taking the same derivatives of the other factors. Thus 



□□ 



□□■ 



-'-'mn-'-'nm 



-'-'tm J-'mn -'-'ni 
J-'UJ-'mmJ-'nn 

etc., 



□ □ 



1 



L 



□□ 



rm 
1 



L r . 



X 



-LJm.n-LJ 



X 



mn ± -"rim 
-t-'mm 



LppL 



U-L-'mm 



(3.60) 



and therefore 










i {2 Hx) = 


( l an ° 


x 


-'-'mn-'-'nm 


(@x -- 


V167T 2 L nn J 


-'-'mm 


-^£mn \ X ) = 


( l an ° 


x 


<-, J- / l.m^ J mn^ J n£ 
^-'(.(.J-'mm 


V16tt 2 L nn J 



Ix = X n ) + subleading, (3.61) 

Ix = Xi) + subleading, (3.62) 



etc., etc. 



The first factor in all these interference terms evaluates to 



1 DD / l 



16tt 2 



r n ' a n 



12 aJK + K, 

v 2 {al + rlf 



(3.63) 



hence 



In the vicinity of the n th instanton Z$ n (x) ~ 2 , ^ lmil 



-^/?™.« \ X J 



D 4 



(3.64) 



Since the standalone instanton's density is In (x) ~ l/o 4 , the interference terms change the 
effective charge radius or rather 



— [due to I^> only] = // d 4 xd 

Pn J J 



ix , I^(x) xI^(x') 4/5 



near X n 



\X — X 



/|2 



(3.65) 
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by 



A n [duetoZ;L] = // 



yjj *Xj Kaj >Xj 



near X n 



A„ [due to Xf mn ] 
etc., etc. 



I / 1 9 

ry ryf _ 

r(l)/^ v r(3) 



a 8 1 

~ — 7 X — - X 



1 



1 



a 4 a 2 L> 2 L> 2 



™> (3-66) 



ax ax -r^- 



X 






l\1 



all a 

~ 75 x 71 x 73 -73- (3-67) 



a 2 " a 4 



L> 4 L> 4 



near X n 



Therefore, the two-instanton interference terms Z^ n affect the Coulomb energy by the amount 
comparable to the direct repulsion (between point-like instantons) but the effect of interference 
terms involving three or more instantons at once is suppressed by extra powers of a 2 /D 2 and 
becomes negligible in the limit of small instantons separated by large distances. Since we have 
already verified similar behavior of the non-abelian energy of the multi-instanton system, this 
completes the proof of our theorem. Quod erat demonstrandum. 

To conclude this section, we would like to derive explicit formulae for the two-instanton 
interference terms (near one of the two instantons) and hence the Coulomb interaction energy 
due to interference. Starting from eq. (3.61) we have the first factor on the RHS evaluated 
in eq. (3.63), so let us evaluate the second factor. For x M = X^, eq. (3.46) gives us 



L mn (x = X T ) = L nm (x = Xj) = (X m - X n )^X 



Yl a £rn a en + 8 ™ fl » X \ ^(vlnVn) 



£^m,n 



(3.68) 

where the first term on the RHS happens to vanish -- cf. eqs. (3.34) -- while the second 
term is 0(a A /D 2 ), which is much smaller than the 0(a 2 ) third term. Thus, 



J-tm.r}. \ X vY n 



J^n.m.\ X vY ri 



a m a n x m^yj 



(3.69) 



while 



hence 



and therefore 
Xi?2 (x near X„ 



^vavayX ^\-n) ' K*-m ^*-n " ! a 



\V Y \ l 



-L-'mn-L-'nm 



[X 



X„ 



a m a n ^ r \}imVn) 



4I„ 



X I 2 



3a^ tr 2 (yt J n j ) a 4 n (a 2 n - r 2 ) a 2 m tr 2 (yly n ) a 2 n - r\ 



(3.70) 



(3.71) 



2 1 V Y" 12 



[ai + r. 



2\5 



n ' nJ 



n\v Y 12 



«n + r n 



xl£\x). (3.72) 
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It is easy to check that such two-instanton interference terms do not affect the net charge of 
the n th instanton — 

J d A xl£l(x) = (3.73) 

near X n 

(2) 

- but they move the charge distribution close to the center X n since Xmn is positive for 
r n < a n but negative for r n > a n . Consequently, they decrease the effective charge radius 
p n , which corresponds to A n > and hence positive extra Coulomb energy. Specifically, 

A, [due to X 2 mn ] = [[ d A x dV -^ {X) X ^ {X>) 



I / 19 



instanton #n 



1 I 4 ^^(^j x ff dr r 3 dr l r >3 _ & n x ««(«« r 'j x 



Q m tr {y m yn) //xJj r y3 a n a n\ a n r 

~r~ /\ I I III I til I ~ ~ ~ r s\ ; ~ r~„ • > , _ _. 

X m -X n \ 2 JJ K + r 2 ) 4 (a 2 +r' 2 ) 5 max(r 2 ,r' 2 ) 

o 

" + sa| x |x„-x,| 2 • (3 ' 74) 

Note however that this A n is a correction of the instanton's charge radius or (rather 1/p 2 ) 
starting from the standalone radius d n instead of the original a n . Consequently, there is an 
additional correction 



4/5 4/5 „ _4/5 _ _4/5 



Ar™ = ^ - ^ * -^ x (Si - aj) = -^ x JJ |ot,|». (3.75) 



The a^ n are spelled out in eqs. (3.35); to the leading order 

VuAXm - X,, 



<n » o§L = "m«n x "" ™_ ^ - x §trfe/ TO (-ir a )) (3.76) 



'-m vv n 



Using rjiZf]J x = <5 a6 cVA, we obtain 



(1) 1 2 ^ Q m a n 



mi;,;,,,!- ~ . - : , x i Y J W{ylyn{-ir a )) = 1 - |tr 2 (^ n ) (3.77) 



'■fimn I I y _ y 1 2 



and hence 

.extra ._ - \ - °m (4 ~ ^{vLVn)) 



1 j- M* MW> + 0(a2/jD 4^ (3>78) 
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Combining the corrections (3.74) and (3.78), we obtain 



A! 



2 ^ 



tT2 (yLyr, 



2 /n4N 



Fi ^^ a 2 I X — X I 2 



+ 0(07,0 



(3.79) 



and therefore 



TV 

C ~ 4AM 



/V4/5 



L^ n 2 ^ IV _ 



f I / m 1 n 

,X m -X n | 2 V H fl n «m 



x (tr 2 (£yj - 2) 



V 



2 /r>4\ 



+ 0{a 2 /D 



I 

(3.80) 



Finally, let us combine the non-abelian and the Coulomb energies of the system and 
re-organize the net energy into one-body, two-body, etc, terms: 



£total = ^£lbody( n ) + I^£ 



2 body , 



m,n) 



+ I Yl £ 3h ° dy (^rn,ri) + 



mj^n 



different 
Lm.n 



(3.81) 



where 



S lbody (n) = N C M (AM 2 x (X*)' 



AM 2 



x at + 



5AM 2 



x 



(3.82) 



g2 body. 



N r 



m,n) 



2AM \X m — X n \ iD 



( A 2 M 4 x a 2 m a 2 n x tr 2 (yly n {-iN mn ■ ?)) \ 

+ 1 - - I -^ + -f > - ^ 2 ^< f 

1 



. ., + -r X (tr 2 (^l/J - 2) 

V a n a mJ 



\ 



+ — 



a 



D 2 XM 2 D 2 



£ 3body (£,m,n) = o( 



Nm 2 



N r 



V AM 2 £> 4 A 2 M 3 L> 4 J ' 



etc., etc. 



J 

(3.83) 

(3.84) 



When comparing magnitudes of the non-abelian and the Coulomb terms here we have used 
a ~ l/(vAM). In fact, we may be more precise by noticing that for small instantons 
distant from each other, the one-body potential energies for the instanton radii o^ are much 
larger than the two-body, etc., interactions between different instantons. Consequently, in 
the minimal-energy or near-minimal-energy configuration of the multi-instanton system, the 
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instanton radii will be close to the equilibrium radius of a stand-alone instanton, 

a n = a + 0(a 3 /D 2 ) (3.85) 

where 

minimizes the one-body potential energy (3.82). Plugging this equilibrium radius into 
eq. (3.83), we obtain a simpler formula for the two-instanton interaction energy: 



2/V 
£ 2hod y(m,n) 



5AM \X m — X n \ 4D 



2 + il2 [yLyn) + tr 2 ^y n (-iiv mn -f)j 

(3.87) 

Note that the expression inside '[•••]' is always positive, so the two-body forces between the 
instantons are always repulsive, regardless of the instantons' 577(2) orientations. However, 
the orientations do affect the strength of the repulsion: two instantons with similar orienta- 
tions repel each other 9 times stronger then the instantons at the same distance from each 
other but whose orientations differ by a 180° rotation (in SO (3) terms) around a suitable 
axis. This fact will be at the core of our analysis of instanton crystals in subsequent sections. 



4 Linear Chains of Instantons 

Consider an infinite set of instantons arranged in a crystalline lattice. In light of eq. (3.87), 
the orientations of instantons at different lattice sites affects the system's energy just as 
much as the lattice geometry. In this section we shall see how this works for a simple one- 
dimensional low-density lattice -- i.e. an infinite chain of small instantons located along a 
straight line at 

X% = (nD, 0,0,0), neZ (4.1) 

where the lattice spacing D is much larger than the instanton radius a. 

Before we go any further, we need a reason why the instantons would form a ID chain 
rather then spread out into all 3 flat space dimensions like holographic baryons (or real-life 
baryons) in the compressed nuclear matter. In principle, we may sidestep this question by 
simply freezing the X% moduli of the instantons and allow only their orientation moduli y n 
to vary, but then it would be hard to place such a system in any kind of a physical context. 
Alternatively, we can make 5D gauge coupling depend on the x<i and x% coordinates as well 
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as the holographic dimension x^ 

8vr 2 



N c \m(\ + M\x\ + M 2 xj + M\x\ + 0(M 4 x A )\ M 4 = M, (4.2) 



g1(x 2 ,x 3 ,x^ 
which creates an effective potential for the instanton centers X n , 



-i\ 



£ (1) (ra) = N c \MMl x (X 4 ) 2 + N C \MM 2 x (X 3 ) 2 + N c \MMl x (X 2 ) 2 + 0(X c AM 5 X 4 

(4.3) 
At low instanton densities, this potential makes the instantons line up along the x 1 axis, 
while at higher densities, the repulsion between the instantons becomes stronger than this 
potential and they spread out into the other dimensions. We shall address such spreading 
out in the following section §5, but for the moment let us focus on the low-density ID 
chains (4.1). 

In a bigger holographic picture, a g$ that depends on three coordinates £2,3,4 indicates 
that the flavor branes are curved in all 3 of these directions, 

[/[the energy scale] « U + ^ + ^ + ^ + 0(X 4 / J R 3 ), (4.4) 

which is quite different from the usual holographic setup where the energy-scale coordinate U 
depends only on the x 4 . But in this article we are not going to explore the geometric aspects 
of such curvature. Instead, we shall simply use the gauge coupling (4.2) in an approximately 
flat space (for distances <C 1/-W) as a tool to put instantons into a ID lattice. In our next 
article [41], we shall make g 2 depend only on x 3 and £4 to put instantons into a 2D lattice 
in the (rci,^) plane, and eventually go back to g§ depending only on the £4 to make 3D 
instanton lattices. 

For simplicity, let us assume that the g§ is much less sensitive to the x 2 and £3 coordinates 
than the £4, thus Mf, M% <^C M\. In this limit, the Mf and Mf parameters give rise to the 
£ 2 and £3 terms in the one-body instanton potential (4.3), but their effect on the instanton 
radius a or the two-body forces between the instantons may be neglected (compared to the 
effect of the M\x\ term). Therefore, the net two-body forces between the instantons remain 
approximately as in eq. (3.87). 

Moreover, for the ID lattice geometry (4.1) of the instanton centers, we have |X m — X n | 2 = 
D 2 x (m — n) 2 while N mn = (±1, 0, 0). Consequently, the net energy of the instanton chain 
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as a function of the instantons' orientations y n may be summarized as 

^ = 5A^ X £ T^hs-2 >< [I + ^(vlVn) + ^{vlVni-ir,))} . (4.5) 

To minimize this energy, each pair of instantons m and n wants to have yj„y n to be a linear 
combination of ir 2 and 2T3; in 5*0(3) terms, this corresponds to a relative rotation (between 
the 2 instantons) through a 180° angle around some axis _L x\. Alas, this cannot be achieved 
for all instanton pairs at once; indeed, if we minimize the energies of the (n, m) and the (n, £) 
pairs, then the energy of the (m, £) pair would be maximal instead of minimal: 

If yly m = m x -T & y\y t = in 2 ■ f for ni,n 2 -1 xi 
Then y} n y e = {n\ ■ f)(n 2 ■ f) = (n\ ■ n 2 ) + i{n\ x n 2 ) ■ f for ni x n 2 || x 1 
=> tr 2 (y^) + tr 2 (y^(-iri)) = 4(n x • n 2 f + A(n 1 x n 2 ){ = 4 (maximum). 

(4.6) 

So our best bet is to first minimize the energies of the most expensive pairs — the nearest 
neighbors (m,m + 1) — and then worry about the less expensive pairs of instantons. Thus, 
we want 

Mm: y\ny m+ i = cosipm x («r 3 ) + sin ip m x (ir 2 ) for some angle ip m , (4.7) 

and the most general solution for this set of equations (modulo a common SU(2) symmetry 
of all the y m ) is 

I ±[cos0 n x 1 + sin0 n x (n~i)] for even n, 

y n = exp(i0 n Ti) x (ir 3 ) n = <j (4.8) 

I ±[cos0 n x (ir 3 ) + sin0 n x (ir 2 )] for odd n, 

for some angles (f> n . Indeed, for any set of the angles ra we have 

ylnVm+i = (-^s) m x exp(i(0 m+1 - m )<Ti) x (ir 3 ) m+1 

= exp((-l)^(0 m+1 - m )n) x (ir 3 ) (4.9) 

= cos(0 m+ i - <f> m ) x (ir 3 ) + (-l) m sin(0 m+1 - m ) x (ir 2 ), 

which agrees with eqs. (4.7) for ip m = (— l) m (0 m +i — m ). Clearly, any set of ip m can be 
obtained for suitable n , so eqs. (4.8) indeed describe all the solutions to the eqs. (4.7). 
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Eqs. (4.8) for various angles <p n define a big family of instanton configurations. Surpris- 
ingly, all these configurations have exactly the same net energies! Indeed, for all sets of 4> n , 
all instanton pairs (m, n) with odd m — n have minimal energies while pairs with even m — n 
have maximal energies: 

for odd m-n, y ] m y n = cos(0 n - <f> m ) x (±ir 3 ) + sin(0 n - <f> m ) x (±ir 2 ) 



2 



l + o + o = \, 



(4.10) 



for even m-n, y\ n y n = cos(0 n - <j> m ) x (±1) + sin(0 n - m ) x (±in] 
=> I + tr 2 (i/^ n ) + tr^yj-in)) = 

= | + 4cos 2 (0 n -0 m ) + 4sin 2 (0„-0 m ) = f . 

Consequently, regardless of the angles <f> m , the net energy per instanton of the ID lattice is 

a /- fixed n 1 

£per instanton = §\M D 2 X Z_/ ( m _ n )2 X 2 + tr l^m^nj + tr (2/mi/n(~ ^ r l) J 

N r sr^ 1 ^ f° r odd 



£ £>' = 



5XMD 2 ^ £ 2 I 9 frir _ VPT1 /; , , „ . 

£=m-n^0 ^2 t0r eVeI1 *-' (4.11 



1^1 9 V 1 

9 Z^ 72 "" O 2-/ 72 



N c / 1 ^ 1 9 

5AMD 2 X I 2^¥ + 2 

\ odd i even ^0 



X - X 



5XMD 2 \2 4 2 12 : 

Now let's consider the instantons' orientations y n from the group-theoretical point of 
view. Each y n is an SU(2) matrix, but its overall sign is irrelevant, so we only care about 
the y n (modulo sign), which belongs to the 5C/(2)/Z 2 = SO(3). In a generic lowest-energy 
configuration (4.8) of the ID chain, the instanton orientations span an SO (2) x Z 2 subgroup 
of the 5*0(3) corresponding to the rotational symmetries of a cylinder — rotations through 
arbitrary angles around the x\ axis, and 180° rotations around axes _L x\. The y n alternate 
between the two types of rotations, but apart from that they generically do not follow any 
regular patterns. 

However, the family (4.8) also contains some regular patterns in which the orientations 
y n (modulo sign) follow a repeating cycle of finite length p; moreover, the values of y n span 
a discrete subgroup of the cylindrical symmetry 5*0(2) x Z 2 . Here are some examples: 
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• The anti-ferromagnetic chain, with 2 alternating instanton orientations: 

Vcvcnn = ±1, Voddn = ±*T 3 . (4-12) 

In this configuration — which obtains for (f) n = -- the y n (modulo sign) span a Z 2 
subgroup of the SO (2) x Z 2 . 

• Period = 4 configuration spanning the Klein group of 180° rotations around the 3 
Cartesian axes 7 : 

Vn=0 (mod 4) = ±1, Un=l (mod 4) = ±^3) Un=2 (mod 4) = ±Tl, Vn=3 (mod 4) = ±^2- 

(4.13) 

• Period = 2k = 6, 8, 10, . . . configurations spanning prismatic groups Z& x Z 2 : 

7rn Tin 

2/evenn = COS — X 1 + Sill — X (iTl), 

( i\ ( i\ ( 4 - 14 ) 

ir(n — 1) ... . 7i(n — 1) v y 

1/oddn = cos — — — x (ir 3 ) + sin — — — x (it 2 ). 

• Period = 2k = 6, 8, 10, . . . configurations spanning dihedral groups D^k-, which obtain 
for (f) n = n x (ir /2k), i. e. 

Tin .Tin 

Vevcnn = COS — - X 1 + Sill— X (in), 

2k 2k u 15) 

Voddn = cos ^T x ( ir s) + sin ^T x ( ir ^)- 

There is a wider class of regular configurations — we shall call them link-periodic — in 
which the y n themselves are not periodic, but the relative rotations y\ l y nJrX between nearest 
neighbors follow a periodic pattern. In terms of the (p n angles, this corresponds to periodic 
differences 4>n+i ~ (pn, for example 

I <z> + 9 for even n, 
<t>n = rnp - \{-l) n 9 =*► 0„ +1 - <f> n = r (4.16) 

I ip — 6 for even n. 

For rational (p/ir this pattern produces a periodic array of instantons' orientations y n — for 
example, the dihedral cycle for tp = tt /2k and 9 = 0, or the prismatic cycle for tp = it /2k and 



7 As a subgroup of SO (3), the Klein group is abelian and isomorphic to Z2 x Z2. But its covering group 
in SU(2) is non-abelian group and isomorphic to the group of unit quaternions ±1, ±i, ±j, ±/c. 



53 



6 = —<p/2. But for irrational </?'s the orientations y n themselves do not have a finite period; 
instead, they wind irrationally around the cylinder group £77(2) x Z 2 . 



All in all, the ID lattice of instantons has a huge degenerate family (4.8) of lowest-energy 
configurations. However, this degeneracy obtains only in the M2,M^ <C M4 limit where the 
5D gauge coupling (4.2) depends mostly on the holographic coordinate X4 and only a little 
bit on the X2 and x%. For the larger M2,M^ ~ M4, the degeneracy is lifted and there is a 
unique lowest-energy configuration of the ID lattice, namely the link-periodic array (4.16) 
whose ip and 9 parameters depends on the M 2 /M 3 ratio. 

To see how this works we note that for M 2 ~ M3 ~ M4, all 3 coordinates ^2,3,4 transverse 
to the ID lattice play similar roles, and the 3 terms M\x\ + Mfx 2 + M^x\ inside '(•••)' 
in eq. (4.2) give rise to similar contributions to the non-abelian energy of a multi-instanton 
system: 

£«* = N c \Mx jd A xI{x) x J2 M l( x n) 2 ( 4 - 17 ) 

(where M\ is the same as M 2 ). Consequently, the two-body interaction energy between the 
instantons becomes more complicated then in eq. (3.42). Instead, generalizing eq. (3.39) we 
obtain 

st a = n c \m J2 K\Y.^ x n) 2 + ¥l) + E^f ( 4 - 18 ) 

/i=2,3,4 L n n=£m ) 

and hence to the leading order in expansion (3.35), 
where 



WaW = V N C \MM 2 (X») 2 + \N c \M{Ml + M 2 + M 2 ) x a 2 (4.20) 



and 



4 2 iKn) = 2| /"^ x |2 x £ iV c AMM 2 x(^iV^trfe„(-^))) 2 . (4.21) 



1 '" * ' " M=2,3,4 



Note that the coefficient of each a 2 in the one-body potential (4.20) depends on the M 2 
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and the M\ as well as the M\ = M 2 . Consequently, when we combine this potential with 
the Coulomb one-body potential for the instanton radius a, the equilibrium radius ends up 
smaller than in eq. (3.86), namely 



2/5 



A 2 M 2 (M| + Ml + Ml) ' 



(4.22) 



Plugging this radius into eq. (4.21) and combining with the Coulomb two-body interaction 
term from eq. (3.80), we obtain the net two-body force in the background with M2 ~ M3 ~ 

M A = M: 



£%*(rn,n) 



2N r 



5AM \X m -X n \ 2 



( \ + tr 2 (yly n ) 

+ J2 C » x (W HvlVni-ir 11 )) 

\ M=2,3,4 



(4.23) 



where 



C„ 



dcf 



Ml 



M 2 + M 2 + M 2 ' 



Ca + Cq + Co 



1. 



(4.24) 



For instantons located along a straight line, this formula can be simplified a bit using 
N v mn = (±1, 0, 0, 0) and hence 



il v K m te(yly n (-ir a )) = < 



±tr fc/nH r1 )) for/x = 4, 
±tr(yt i |/ n (-ir 2 )) for// = 3, 
±tr(^ n (-ir 3 )) for// = 2. 



(4.25) 



Consequently, the net interaction energy of the ID lattice is 

N, 



cuit 



^net 



2 ^ 



(2)( 



x ^ Q(m,n) 

X 



E 



5AML> 2 ^ (m- n) 2 



(4.26) 



where 



Q(m,n) = \ + tr 2 (yiy n ) + Ca tr 2 (yly^-tr 1 )) (4.27) 

+ C 3 tr 2 (yly n (-tr 2 )) + C 2 tr 2 (yly n (-ir 3 )) . 



Without loss of generality we assume Ma > M3 > M 2 and hence Ca > C3 > C 2 - - oth- 
erwise we simply re-label the coordinates #2,3,4. Thus, the lowest-energy relative orientation 
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of an instanton pair (m, n) is yl n y m = ±ir 3 . If this orientation is unachievable, the next best 
bet is a linear combination of ir 3 and ir 2 . If that is also unachievable, the third-best choice 
would be a linear combination of all three ir 1 ' 2,3 ; in SO (3) terms, this corresponds to the 
180° rotation around a generic axis. Finally, rotations through other angles would be the 
least attractive option which the instantons #m and #n would rather avoid — unless they 
are forced to it by interactions with the other instantons between m and n. 

In this setting, it is not intuitively obvious how to balance the energy cost of nearest- 
neighbor interactions versus next-to-nearest neighbors and more distant instanton pairs. 
Instead of intuition, we have performed a computer experiment using a lattice of 200 SU(2) 
matrices y n . In each run, we start with the y n being independent random elements of the 
SU(2) group, and then let them evolve towards a minimum of the energy function (4.26) via 
the relaxation method. That is, we let the y n evolve with time according to 



dy n (t) S£, 



hit 
net 



K x — ^ (4.28) 

dt 5y n v ; 

where K is a constant mobility factor and the derivative with respect to an SU(2) matrix 
y n is defined as 

^- = yn((-ir)-V s £{y n ^y n (l + is-f))\. =o y (4.29) 

Each run ends when all the y n seem to converge to an equilibrium configuration and their 
derivatives (4.28) become very small. We made many runs for different C3 and C4 parameters, 
and here is what we have found: 

• In all equilibrium configurations for backgrounds with C 3 < C 4 , all twists yjjy n+ i 
between nearest neighbors have form 

ylVn+i = cos^n x ir 3 + sin^„ x ir 2 + (tiny) x ir 1 + (tiny) x 1, (4.30) 

where the tiny coefficients of 1 and ir 1 are artefacts of imperfect convergence to equi- 
librium (they get smaller when we allow more time for convergence). 



o 



The backgrounds with M 3 = M 4 =^> C3 = C4 are trickier due to symmetry between 
ir 1 and ir 2 components of the 577(2) matrices. (Or all three ir 1,2 ' 3 matrices for M 2 = 
M3 = M4 =^> C2 = C3 = C4.) In equilibrium configurations for such symmetric 
backgrounds, the nearest-neighbor twists y„y n+ i generally do have ir 1 components. 
However, all such components can be eliminated by a suitable global symmetry of the 
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system, y n — > y n x G, same G G SU(2) for all y n , and then all the nearest-neighbor 
twists y^y n+1 take the form (4.30). 

For most nearest- neighbor pairs, the angles ip n in eqs. (4.30) take values ±<£> (mod 7r) 
where y2 depends on the 6*2,3,4 parameters. Moreover, the ± sign here tends to alternate 
between odd and even n. However, some pairs do not follow these rules. 

In general, our runs end up with patterns of (—l) n if) n that look like this: 

(-l)^n 



-*■ n 



(4.31) 
Physically, this corresponds to the ID lattice having two degenerate ground states, 
one with ip n = + (— l) ra y? and the other with ip n = — (— 1) T V- in our numerical runs, 
some parts of the lattice converge to one of these states while other parts converge to 
the other state; altogether, we end up with several domains separated by walls. In a 
perfect simulation, the walls would eventually move towards each other and annihilate, 
so the whole lattice would end up in the same ground state. But this process is so slow 
that our numeric simulation stops before it barely begins, thus we always end up with 
multiple domains instead of a single ground state for the whole lattice. 

* The bottom line of our numerical simulations is that the ground state of the ID in- 
stanton lattice always has a link-periodic instanton orientations (4.16) (up to a global 
symmetry, if any) for a periodicity angle tp that depends on ratios of parameters 
M2 : M3 : M4. Curiously, it actually depends only on the M2 : M3 ratio regardless of 
the M 4 (as long as M 4 > M 3 > M 2 ). 

Once we know the ground state of the instanton orientations in a ID lattice is link- 
periodic, we may calculate the periodicity angle <p analytically. Indeed, applying eqs. (4.10) 
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to the link-periodic array (4.16) of instanton orientations, we have: 



For even n — m, y m y n 



cos((n — m)ip) x i n m + sin((n — m)(p) x i 



n+m+l 1 



Q(m,n) = \ + 4cos 2 ((n — m)(p) + AC4 sin 2 ((n — m)(p) 



For odd n - m, y] n y n 



Q(m,n) 



2 

§ + 2C 4 + 2(1 -C 4 = C 3 + C 2 ) xcos(2(n-m)v9), 

cos((n - m)y? ± 0) x i"- m r 3 + sin((n - m)y? ±0) x i n+m r 2 
\ + 4C 2 cos 2 ((n-m)ip± 6) + 4C 3 sin 2 ((n - m)ip ± 0) 
\ + 2(C 3 + C 2 ) - 2(C 3 -C 2 ) xcos(2(n-m)v?±2fl), 



(4.32) 



(4.33) 



where the ± sign is (— l) n . To calculate average interaction energy per instanton in the 
lattice, we should average Q(n, m) between odd and even n for fixed £ = n — m; indeed 



c interaction 
per instanton 



N r 



5XMD 2 



x 



fixed n 



^ Q(m,n) ^ BV<Taft " 

f J In — mn\2 



m=£n 



(n — my 



A', 



5XMD 2 



£^0 



QW 



Thus, 



(4.34) 



for even £, Q(£) 



+ 2C 4 + 2(C 3 + C 2 )cos(2^>), 



for odd £, Q(£) = \ + 2(C 3 + C 2 ) - 2(C 3 -C 2 ) xcos(20) xcos(2^>), 
and therefore 



even <?^0 



<7 interaction c 

C pcr instanton ~ tj\lU"n2 



12 



+ 2(C 3 + C 2 ) x ]T 



(4.35) 



Vcvcn £^0 



cos(2£ip) IT 2 . . 

^—^ = 7tU + 111? 

£ 2 12 m * 



+ (| + 2C 3 + 2C 2 )x(^i = ^ 



vodd^ 



2(C 3 - C 2 ) cos(20) x J2 



vodd^ 



cos(2£y9) 

£2 



TV' 



-j-*m 



\ ((where the sums are evaluated assuming \ip\ < | )) J 
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/ 



N c 
5AML> 2 



\ 



yx(l+C 3 + C 2 -(C 3 - C 2 ) cos(2fl) 
2n\p\x(c 3 + C 2 - (C 3 - C 2 ) cos(29) 

+ Vx(c 3 + C 2 ) 

Minimizing this expression WRT the </? and produces 4 degenerate minima, namely 

C 2 



(4.36) 



/ 



(1) 
(2) 
(3) 
(4) 






7T 

+ 2 X 

7T 
X 

2 

7T 
X 

2 



7T 

+ 2 X C 2 + C 3 



C 2 + C*3 

c 2 


C*2 + C,3 

c 3 


C 2 + C*3 

c 3 












— 
2< 

7T 

2' 



(4.37) 



however the last two minima are physically equivalent to the first two. Thus, in agreement 
with our computer 'experiments', the ID instanton lattice has 2 degenerate ground states 
related by the <p — > —(p symmetry. The value of \ip\ also agrees with our 'experimental data': 







ix M 2 

0, ip = ±— x — k ^ 

' r 2 M 2 + M 2 



regardless of the M 2 



(4.38) 



(as long as M 4 > M 3 > M 2 ). In particular, for rational ratios of M%/M 2 , the angle p> is 
rational (in units of 7r) and the instanton orientations y n become periodic. Specifically, for 



M 2 P ,, x 

— 2 = - , p < q, gcd(p, q) 
ivi 3 q 



(4.39) 
dihedral sub- 



the y n (modulo signs) repeat with period 2(q + p) while spanning the D 2 ( q+P ) 
group of the 5*0(3) = S r £/(2)/Z 2 , for example see eq. (4.15) for p — 1 and q — k — 1. 

Two particularly interesting Mf : Mf ratios need special handling, Mf = Mf and Mf <C 
Mf . For M 2 = M 3 the background has a rotational symmetry in the x 2,3 plane; for the ID 
instanton lattice, this translates into the U(l) symmetry between the ir 2 and ir 3 directions in 
the SU(2). Consequently, instead of two discrete ground states (4.38) we have a continuous 
family: 

71 

ip — — , 9 = anything. (4.40) 



For all configurations in this family, the instanton orientations y n (modulo signs) repeat with 
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period 4 while spanning the Klein groups Z 2 x Z 2 of 180° rotations around 3 mutually _L 
axes. For 9 = — f the y n are spelled out in eq. (4.13); for other values of 9 we have similar 
cycles related by the U(l) C SU(2) symmetry. 

Finally, for the very asymmetric background with Mf <C M 2 , the two ground states 
(4.38) become indistinguishable as <p — >■ ±0. In this limit, the instanton orientations form 
the anti- ferromagnetic order (4.12). 

5 Instanton Zigzags 

In this section we consider the first step in the transition between ID and 2D instanton 
lattices — a lattice with one infinite dimension x 1 while the other dimension x 2 has just two 
layers. Since the instantons repel each other, the two layers should be staggered in the x 1 
direction, so the whole lattice looks like a single zigzag-shaped chain: 

x 2 



*x SK ^ 



.vr v^ 



•^ £ 



.vr vr. 



x 



(5.1) 

In such a zigzag, the instanton j^n has its center at 

X\ = n x D, X 2 n = (-1)" x e, X 3 n = X n 4 = 0, (5.2) 

for some parameters D and e; we shall refer to D as the lattice spacing and to e as the zigzag 
amplitude. 

To make the instantons form a zigzag lattice, we use the same background 



8vr 2 



g 2 (x 2 ,x 3 ,x 4 ) 



N c XM(l + M 2 x\ + M 2 x 2 3 + Mfa^ + 0(MV)Y (4.2) 



we have used in the last section to make a straight chain but increase the instanton density 
p = 1/D until the repulsion between the instantons pushes them away from their neighbors 
in some transverse direction. For M 2 < M 3 < M 4 , the initial breakout from the straight 
line is going to be in the ±x 2 directions, with opposite signs for the nearest neighbors (to 
increase the distance between them), hence the zigzag geometry (5.2). For higher densities 
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- and hence stronger repulsion between the instantons - - we will get more complicated 
lattices: multiple layers in the x 2 direction, and eventually breaking out into the x 3 and x 4 
dimensions. However, the first geometry right after the straight ID lattice is going to be the 
zigzag. 

Since our formulae for the instanton interactions presume distances \X m — X n \ much 
larger than the instanton radius a, we need the transition from a straight ID chain to a 
zigzag to happen at a critical lattice spacing D c ^> a, which calls for 

M 2 <t: M 4 = M. (5.3) 

Indeed, consider an instanton that have moved a little bit away from the x 1 axis in the x 2 
direction. The restoring force on this instanton due to x 2 dependence of the gauge coupling 

is 

F rest = -X 2 x 2N C \MM 2 (5.4) 

while the repulsion from the neighboring instantons pushes it further away with the net force 

t? d n( N * 1 ^ X 2 ( N c \ 

F «* = ~dX° [AM X 0(D 2 )+X 2 ) = + D~ 4 X ° \XM ) ■ (5 ' 5) 



At the critical lattice spacing, the two forces balance each other, hence 

0(1) 
\ 2 M 2 M 2 



g 4 _ Q(l) (56) 



At the same time 



2/5 0(1) 



a X 2 M 2 {M 2 + M 2 + M 2 ) \ 2 M 2 M 2 ' (5 ' 7) 

so to assure D c ^> a we need M 2 <^ M4 = M. 

Note however that while we need M 2 ^C M4, there are no restrictions on the M3 parameter 
except M 2 < M3 < M4. Consequently, we may consider zigzags in a variety of backgrounds 
with different M 3 /M 4 ratios - - ranging from M 2 = M 3 < M 4 to M 2 < M 3 = M 4 - 
and hence different ways in which the forces between two instantons depend on their relative 
orientations. For small M3 <^C M4, the two-instanton interactions are spelled out in eq. (3.87) 
while for finite M3/M4 ratios we need to use a more complicated formula (4.23) (with C 2 ~ 
but finite C3 and C4). For a zigzag where all the instantons lie in the x 1,2 plane and hence 
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all the N^ n have form N^ n = (*, *, 0, 0), we may simplify this formula using the identity 
E (^mntr(^ n (-zr a ))) 2 = J2 ^ 2 {yly n {-iT a )) (forAC = (*,*,0,0)). (5.8) 



/i=2,3 a=l,2 

Consequently, eq. (4.23) can be rearranged as 



Zigzag] = 7TT7 2. |y_y.|2 ^ 9a ) 



5AM *—f \X m — X n , 



where Q z (m,n) = \ + tr 2 (^ n ) + C 3 ^ tr 2 (^„(-ir a )) 

a=l,2 

+ (1 - 2C 3 ) x tr 2 {yly n (-ir ■ N mn j) . (5.9b) 

Note that in general, the force between two instantons is not central — it depend not only on 
the distance between the instantons and their relative orientation y^y„ of their isospins, but 
also on the direction N mn of their separation X n — X n in space. However, in the background 
with M 3 = M4 — which is the background we have used in our previous paper [1] (albeit in 
different notations) -- C3 = |, the last term in eq. (5.9b) goes away, and the force becomes 
central (but isospin-dependent). 

The forces between instantons will be important in our next article [41] where we shall 
allow them to move in the x 1,2 plane seeking the lowest-energy lattice. For the moment, we 
simply assume that the instantons somehow form a zigzag of some particular lattice spacing 
D and amplitude e and focus on their orientation moduli y n . Our immediate goal is to find 
the lowest-energy configuration of the orientations y n as a function of the e/D and M3/M4 
ratios. 

A priori, we do not know if the orientations — or rather the relative orientations y^y n+ i of 
the neighbor instantons — follow any particular pattern or patterns. To find such pattern(s) 
we used numerical simulations: starting with completely random y n , we let them evolve 
towards a minimum of the energy function (5.9) according to eqs. (4.28); when the evolution 
seems to stop, we look at the relative orientations of the nearest neighbors yljjn+i to see 
if they form any recognizable pattern. Repeating this procedure for different combination 
of the e/D and M 3 /M 4 parameters, we saw five distinct patterns of orientations (or rather 
four distinct patterns and one confused mess). Figure 5 diagrams the distribution of these 
patterns in the parameter space. 



62 



»••••' 



• ••••I 



M3/M4 

M 

1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



0.0 



'•••••'•••••' 



•••••••••• 



;•••• •!•••••! 






• •••• 

• ••••! 



•••••••••• 



•••••••••• 



•••••••••• 



■ J- ■ 



• ••• alaaaaa! 



I 



i* •••••••• 



a aa alaaaaa! 



laa a ••!••• I 



•I* 



•••••••••• 



•••••••••• 



• ••••'• ••••••••••••• • 






• ••••!• ••••!• ••••!•••• • 



I 



-I- 



I 
■ 4- 



• • •••!•••••!• •••••••• •> 

••!•• ••!•• ••!•• 

III! !!•! ■>!>■ 



• • • ••!•••••!• ••••!•••• •! 



-I- 



. 4- 
I 
I 
I 



• ••• 

• ••• 



• ••• 

• ••• 

• ••• 



la • • • 
l??-"-5' 



-^ e/D 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Figure 5: Rough distribution of instanton orientation patterns in the zigzag parameter space. 
This diagram is obtained from a numerical relaxation method. 

Here is the color key to this figure: 

• The red dots on figure 5 denote the anti-ferromagnetic pattern (AF) of instanton 
orientations in which the nearest neighbors always differ by a 180° rotation around the 
third axis, 



±1 for even n, 
±ir3 for odd n, 



same y\ l y n+l = ir 3 for all n. 



(5.10) 



The yellow dots denote another abelian pattern ( ) in which all nearest neighbors 
differ by the same U(l) C SU(2) rotation, but now the rotation angle is < 180°, 



same y] l y n+1 = exp(|0r 3 ) for all n, < 4> < ft- 



(5-11) 
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The blue dots denote a non-abelian link-periodic pattern (NAl) in which the relative 
rotation between nearest neighbors is always through a 180° angle, but the direction 
of rotation alternates between two different axes in the (12) plane, one axis for the 
odd-numbered instantons and the other for the even-numbered. In SU(2) terms, 



v\kVvk+i = exp(fn e -f) 
vlk+iV2k+2 = exp(f n -f) 

for some A, B ^ {A 2 + B 2 = 1). 



in e ■ t = 


= +iAri + iBT2 , 


in ■ t - 


= +iAri — iBr 2 , 



(5.12) 



The green dots denote another non-abelian link-periodic pattern (NA2). Again, the 
relative rotation between nearest neighbors is always through a 180° angle, but the 
direction of rotation alternates between two different axes. However, this time the two 
axes no longer lie within the (12) plane, thus 

ylkV2k+i = iAr i + iBr 2 + iCr 3 , 

(5.13) 

y2k+iV2k+2 = %Ar \ - iBt 2 - iCr 3 , 

where A, B, C all ^ (A 2 + B 2 + C 2 = 1). 

• Finally, the purple dots denote a confused non-abelian pattern NAX in which the near- 
est neighbors seem to differ by rotations through angles <fi < 180° around 2 alternating 
axes. In SU(2) terms, we see something like 

ylkV2k+i = + iAr i + i Br 2 + iCr 3 + D, 

(5.14) 

y2k+lV2k+2 = -lATl + lBT 2 + lCT 3 + D, 

where A, B, C, D all ^ (A 2 + B 2 + C 2 + D 2 = 1), but this picture is rather noisy: 
the coefficients A,B,C,D fluctuate from n to n and from run to run much stronger 
than their analogues in the other 4 patterns (5.10)-(5.13), and even the relative signs 
between A, B, C, D flip every few lattice sites — the long-range correlations are quite 
poor. 

Physically, the disorder in the purple part of the phase diagram 5 indicates a nearby first- 
order phase transition between incompatible orientation patterns. Near such a transition, 
our relaxation method fails to converge to a homogeneous phase that happens to have lowest 
energy for the {e/D, M 3 /M^) parameters in question. Instead, it leads to a mishmash of 
short domains of two nearly-degenerate phases; in fact, the domains are only a few lattcie 
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sites long, so they are significantly distorted by the boundary effects. The resulting mess is 
not even a local minimum of the zigzag's energy function, but once our simulation reaches 
this point, further evolution of the instanton orientations y n becomes very slow. Given very 
long time (and a better numeric algorithm) we would eventually see the domains of each 
phase growing longer and eventually merging into a single phase -- but in our calculations 
we have simply run out of time and patience before this can happen. 

The bottom line of our numeric analysis is that the instanton zigzag has 4 or 5 distinct 
orientation phases — we can not tell if NAX is a real phase or an artefact of mixed domains 
of other phases that we could not resolve. To answer that question - - and also to get a 
better map of phase boundaries — we need to go back to analytic work. 

Our starting point is the qualitative result of our numerical analysis: all 5 patterns 
(5.10)-(5.14) we have seen — or might have seen — can be fit into a single anzatz: 

y n = exp(mfr 2 ) x exp(z(f + (-l)*f )n), (5.15) 

VnUn+i = cos | cos/3 x 1 — (— l) n cos| sin/3 x %T\ 

+ sin | cos a x ir 2 + sin | sin a x ir 3 (5.16) 

for some angles <fi, a, /3. Indeed, the AF phase obtains for </> = n, a = \\ the phase 

obtains for < </> < ir, a = | , /3 = 0; the NAl phase obtains for < < ir, a = 0, (3 = |; 
the NA2 phase obtains for < < 7r, < a < ~, — ~; and the NAX phase — if it exists 
at all — obtains for generic values of all three angles, < (ft < it, < a < f, 0</3< : |. 

Given the ansatz (5.15), we may calculate the zigzag's net energy as an analytic func- 
tion of the angles 0, a,/3 and parameters (e/D, M3/M4). The calculation is presented in 
Appendix A, and the result is 

CTnLnton = ^^2 x H<f>, «, fr Cs, e/ D) + N C \MMJ x e\ (5.17) 

3 



for T = - + (-1 + (1 -C 3 )(cos 2 acos 2 /3 + sin 2 asin 2 /3)) x S o (0) 

+ (- + cos 2 a - sin 2 /3) x Si(e/D) + (2 - cos 2 a - sin 2 /3) x S 2 (0, e/D) 

+ (1 - 2C 3 )(sin 2 /3 - cos 2 a) x S 3 (e/ J D) 
+ (1 - 2C 3 )(sin 2 /3 + cos 2 a) x S 4 (0, e/D) 
- 2(l-2C 3 )cosasin/3x E 5 {(f>,e/D), 

{5.11 
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where 



£„ = ^LfW, (5.19a) 



« 2 



= -(^^ (5lgb) 

_ siiih((7r-2^)e/£>) 
2 " (7re/D)cosh(7re/D)' {b ^ C) 

^ = crf(W (5 ' 19d) 

cosh(20e/D) _ 20 cosh((7r - 20)e/D) 

4 ' cosh 2 (7re/Z>) vr cosh^e/D) l ' 6J 

sinh(20e/D) 20 sinh((7r - 20)e/D) 

5 '^ cosh 2 (7re/D) vr cosh^e/D) ' { ° * > 

In Appendix A we also minimize J 7 - - and hence the net energy per instanton — with 
respect to the orientation parameters a, (3, and for fixed C% and e/D. In particular, we 
show that the minimum never lies at generic values of the angles 0, a, fi - - which would 
correspond to the NAX phase. This means that the NAX phase does not really exist — its 
appearance in the purple region on the figure 5 is an artefact of poor convergence near a 
first-order phase transition (see below). 

On the other hand, for appropriate values of the e/D and M3/M4 parameters, the zigzag 
energy (5.17) has global minima corresponding to the other 4 phases (5.10)-(5.13) of instan- 
ton orientations. Figure 6 on the next page shows their distributions in the parameter space. 
Curiously, one of the non-abelian phases — the NA1 -- appears in two different regions of 
the parameter space -- one at low M^/M^ and small zigzag amplitudes, and the other at 
higher M3/M4 and higher amplitudes — while the intermediate region is occupied by the dif- 
ferent non-abelian phase NA2. Both phases transitions between the two non-abelian phases 
are second-order. Likewise, the transition between the non-abelian NA2 and the abelian 
anti-ferromagnetic phase AF is also second-order. 

On the other hand, the transitions between the other abelian phase and the non- 

abelian phases NAl and KA2 are first-order. In fact, these are precisely the first-order 
transitions that confused our numeric analysis in the purple region of the figure 5 into 
appearance of the fifth orientation phase NAX. 

Finally, the triple point at (M3/M4) = 1, (e/D) w 0.38 between the AF, NA2, and 
phases is critical. At that point, the transitions between all 3 phases are second-order. 
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M 3 /M 4 




.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 6: Phase diagram of the instanton orientation patterns in the zigzag parameter 
space. This diagram obtains from the mostly-analytical calculation in Appendix A. The 
lines separating different phases indicate the order of the phase transition: a white line for 
a second-order transition and a black line for a first-order transition. 

Indeed, in our previous paper [1] - - where the analysis was limited to zigzags with M 3 = M4 
(in present notations) --we saw a second-order transition between the two abelian phases 
AF and 



Thus far, we have focused on the instanton orientations y n for a given zigzag geometry. 
That is, we have not only assumed that the instanton centers form a zigzag, but we have 
also treated the zigzag amplitude e as an independent input parameter, just like the lattice 
spacing D or the M 2 , M 3 , M 4 parameters of the 5D gauge theory. Indeed, figure 6 shows the 
transitions between different orientation phases as function of independent, freely- adjustable 
parameters e/D and M3/M4. But physically, the zigzag amplitude e is a dynamical modulus 
- in fact, all instanton center coordinates X^ are dynamical moduli — whose value follows 
from minimizing the net energy of the multi-instanton system. 
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So our next task is to find the zigzag amplitude e as a function of the instanton density 
p = 1/D by minimizing the net energy (5.17) of the zigzag with respect to e (in addition 
to the orientation moduli a,/3,<f>). We begin by changing variables from e to £ = e/D and 
rewriting eq. (5.17) as 



<?zigzag 
per instanton 



ir 2 N c p* 
20AM 



x 



.F m (£;M 3 /M 4 ) + 



20A 2 M 2 M 2 
7r 2 p 4 



xf 



where 



J- m (£;M 3 /M 4 



min J-(a, (3, 0; (f = e/D), M 3 /M 4 ) 

a,p,<P 



(5.20) 



(5.21) 



for J 7 from eq. (5.18). Note that although we know J 7 as an analytic function of all its 
arguments, and we know how to minimize it analytically with respect to the a and (5 angles 
(cf. Appendix A), the minimization with respect to the angle has to be done numerically, 
so we do not have an explicit formula for the J- m (£\ M3/M4). Nevertheless, it is fairly easy 
to calculate numerically both the T m function itself and its derivative dJ- m /d^. 

Given eq. (5.20) for the zigzag energy, the minimum condition dE/de = becomes 

equivalent to 

BE 8T„ 20A 2 M 2 M| 

2£ x - r - r ^ = 0, (5.22) 



■X 



7T 2 p 4 



<9£ " d£ 

which gives us the implicit solution - - instead of the energy-minimizing e or £ = e/D as 
a function of the instanton density p, we obtain the density as a function of the energy- 
minimizing £, 

p\£) = 22 (AMM 2 ) 2 x ZL ^ . (5.23) 



K 



For example, evaluating the RHS of this formula numerically for the M3 = M4 background 
- which is equivalent to the model studied in out previous paper [1] - - we obtain 
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Naively, given such a graph, we may simply swap its axes to obtain the explicit solution for 
the zigzag relative amplitude £ = e/D as a function of the instanton density p: 
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(5.25) 
Note however that this graph has two major limitations: First, it spans all the £ > zeroes 
of the derivative dS/dC, = 0, including all the local minima of the zigzag's energy function, 
and even the local maxima, if any; consequently, for some p we have multiple solutions for 
the zigzag relative amplitude £. Second, while going from eq. (5.22) to eq. (5.23), we have 
missed the straight-chain solution £ = 0, which is always a local extremum of the energy 
function. 

Therefore, the actual procedure for finding the equilibrium zigzag shape parameter £ = 
e/D as a function of the instanton density p = 1/D involves an extra step: For each p, we 
list all the zeros of the dS/dC, derivative, including all the £ > values from the graph (5.25) 
and also £ = 0, then for all these zeroes we calculate the energy function £(£) to see which 
zero is the global minimum of energy. Thus, for zigzags in the (M3/M4) = 1 background we 
ultimately find 
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Note two distinct phase transitions as we increase the instanton density p. First there 
is a second-order phase transition from the straight-chain to a small-amplitude zigzag at 
p c i ~ 1. 053V AMM 2 , and then there is a first-order transition to a larger-amplitude zigzag 



phase at p c2 ~ 1.26\AMM 2 . Besides the discontinuous £(p) function, the second transition 
also changes the instanton orientation pattern from the anti-ferromagnetic AF to the abelian 
with < 7T. On the graph (5.26), these orientations patterns are indicated by the colors 
of the £(p) line — red for the AF and yellow for the 

Besides the M3 = M4 background from our previous paper, we have performed similar 
calculations for 100 new backgrounds with M 3 /M 4 ratios ranging from 0.00 to 0.99. We 
have performed similar calculation for zigzags in backgrounds with 100 different M^jM^ 
ratios ranging from to 1. We found that usually there are multiple phase transitions as 
the instanton density increases, but the specific transition sequence depends on the M3/M4 
ratio! Also, for most M3/M4 ratios, at least one of the zigzag phases is non-abelian! For 
example, for (M3/M4) = 0.8 there are 3 rather than two phase transitions: First, there is 
a second order transition from the straight chain to the AF phase of the zigzag, then there 
is a first-order transition to the non-abelian NA2 phase of the zigzag, and then another 
first-order transitions to the other non-abelian phase NA1. Graphically, 
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For another example, for (M3/M4) = 0.2 there are only two phase transition, but this time 
the first transition — from the straight chain to the NAl phase of the zigzag — is first-order, 
while the second transition — to the NA2 phase — is second order. Graphically, 
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Figure 7: Diagram of the zigzag and straight-chain phases in the p 2 vs. M 3 /M 3 plane. The 
straight-chain phase is colored pink, while 4 other colors — red, yellow, blue, and green - 
denote zigzag phases with different instanton orientation patterns. The first-order transition 
between phases are indicated by black lines, the second-order transitions by white lines. 

Altogether, we have found 7 different transition sequences for different M3/M4 ratios! 
But instead of spelling out all these sequences in gory detail, let us simply present in figure 7 
the phase diagram of all the zigzag and straight-chain phases in the p 2 versus M3/M4 plane. 
Here are a few particularly noteworthy features of this diagram: 

• The very first transition from the straight chain to a zigzag could be either first-order 
or second-order, depending on the M3/M4 ratio: for (M3/M4) < 0.63 the transition is 
first-order while for (M 3 /M 4 ) > 0.63 it's second-order. 

• The non-abelian phases NAl and NA2 of the zigzag cover much larger areas of the 
phase diagram 7 than the abelian phases AF and . In particular, at larger densities 
p -- and hence larger zigzag amplitudes -- the instanton orientations usually prefer 
non-abelian patterns. Only the backgrounds with M 3 « M 4 -- such as the model we 
have analyzed in [1] — favor the abelian orientations. 
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• In this section, all the backgrounds have M 2 <C M 3 , so the straight chain phase on 
figure 7 always has the anti-ferromagnetic order of the instantons' orientations. This 
explains the order of the phase transition from the straight chain to the zigzag: For 
(M3/M4) > 0.63 the zigzag phase right after the transition has the same AF orientation 
pattern as the straight chain, so the transition between them is second-order. But for 
(M3/M4) < 0.63 the zigzag has a different orientation pattern — the non-abelian NAl 
or NA2 — so the transition is first-order. 

• At (M3/M4) = 0.885 and (p 2 /XMM 2 ) = 1.56, all 4 zigzag phases seem to meet each 
other at a quadruple point. We do not know if this an un-resolved pair of triple points 
or a genuine quadruple point; a more accurate numerical analysis should resolve this 
issue. 

• At (M3/M4) = 0.69 and (p 2 /XMM 2 ) = 1.75 there is a critical point on the transition 
line between the two non-abelian phases NA2 and NAl: the discontinuous jump of e 
as a function of p disappears, and the transition's order changes from first to second. 

The phase diagram on figure 7 completes out analysis of the instanton zigzag. But the 
zigzag is only the first step in the transition sequence from a ID to a 2D instanton lattice. 
The subsequent steps — from a zigzag to a 3-layer lattice, to 4-layer lattice, etc. -- will be 
explored in our next paper [41]. 

6 Summary and outlook 

In [1] we have investigated the crystalline structure of the holographic nuclear matter. The 
analysis done in that paper revealed the "popcorn transition" of a one- dimensional infinite 
chain of nucleons into a zigzag structure. We have also found out that the phases between 
adjacent nucleons are abelian. In this paper we continued the exploration of ID instanton 
systems and their popcorn transition to the zigzag, but by generalizing the background to 
allow for M 3 7^ M 4 , we found a much richer phase structure, including two distinct non- 
abelian phases. 

In §3 we proved that for low-density multi-instanton systems --in which the distances 
between the instantons are much larger that their radii - - the interactions between the 
baryons are dominated by the two-body forces, while the 3-body forces, etc., are suppressed 
by powers of (radius/distance) 2 . The two-body-force approximation makes for much easier 
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calculations, so in §4-5 we could put the instanton chains and zigzags into richer backgrounds 
than we had in our previous paper. In particular, we could use a less symmetric external 
potential for the baryons, by making the 5D gauge coupling vary as 

-^— = N C XM x (l + M\ xx\ + Ml xx 2 3 + Ml x x\ 0(M 4 ^ 34 )) (1.2) 

with three independent curvature parameters M 2 , M 3 , M 4 . 

For the almost-holographic background M 2 , M 3 <C M 4 , we found that the ID lattice of 
instantons has a highly degenerate family of ground states with many different patterns of 
instantons' orientations. Among them are periodic patterns in which instanton orientations 
span finite subgroups of the S'L r (2)/Z 2 : the anti-ferromagnetic chains spanning the Z 2 sub- 
group (two alternating orientations related by a 180° twist), period = 4 chains spanning the 
Klein group Z 2 x Z 2 , and period = Ik chains spanning the prismatic groups Z& x Z 2 and 
the dihedral groups D 2 k- There are also many link-periodic patterns in which the relative 
orientations y„y„ + i of neighboring instantons are periodic but the orientations y n themselves 
are not periodic. But the vast majority of the degenerate ground states are not periodic 
at all. Raising the M 2 and M 3 parameters lifts the degeneracy. For M 2 <C M 3 ~ M 4 , the 
ground state is the anti-ferromagnetic chain, while for M 2 ~ M 3 ~ M 4 the ground state 
is periodic with a longer period or link-periodic, depending on the M 2 /M 3 ratio. Specifi- 
cally, for M 2 = M 3 the instanton's orientation in the ground state span the Klein group, 
for other rational Mf / 'Ml they span a dihedral group, while the irrational M^/Ml favor the 
link-periodic patterns. 

For the zigzag-shaped chains of instantons we need M 2 <C M 3 ~ M4 to make sure the 
transition from a straight chain to a zigzag happens for lattice spacing 3> instanton radius, 
but we may vary the M3/M4 ratio. Consequently, we found 5 distinct phases shown in fig- 
ure 7: The anti-ferromagnetic straight chain, the anti-ferromagnetic zigzag, the abelian link- 
periodic zigzag (period = 1), and two different non-abelian link-periodic phases (period = 2). 
Figure 7 also shows a rather complicated phase structure: As we increase the instanton den- 
sity, the zigzag amplitude and the pattern of orientations go through a sequence of 1, 2, or 
3 phase transitions, but the number of transition and their orders depend on the M3/M4 
ratio. In particular, the very first transition from a straight chain to a zigzag is first-order 
for M 3 < 0.63M 4 but second order for M 3 > 0.63M 4 . The subsequent transition(s) may also 
be first-order or second-order, depending on the M 3 /M 4: . 

This paper is part of a program focused on exploring the holographic nuclear matter, 
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especially its low-temperature high-pressure phases. Thus far we have worked out the ID 
instanton lattices and their zigzag deformations as a toy model for the transitions between 
3D and 4D lattices. But many interesting questions remain open, and we would like to 
address them in our future papers; here is a short list: 

• The very next step in our program is to follow up the "popcorn transitions" from a ID 
instanton chain into the second dimension beyond the zigzag phases through multiple 
lines of instantons towards a thick 2D lattice. After that, we should set M 2 = and 
study the infinite 2D lattices and their popcorn transitions into into the third dimension 
for M 3 <C M 4 . The techniques developed in this paper -- using the two-body-forces 
approximation for a <C D - - should work all the way up to infinite 3D lattices for 
M 2 = M 3 = 0. Alas, the ultimate popcorn transition from the 3D instanton crystal 
into the fourth space dimension probably happens for D ~ a, so the two-body-forces 
approximation would not be accurate. But we ought to to try it anyhow, just to get a 
qualitative picture of the transition. 

• Our program (starting with [1]) is based on the conjecture that the popcorn transition 
from a 3D instanton crystal to a 4D crystal is holographically dual to the transition of 
the large- A/" c nuclear matter from the baryonic phase to the "quarkyonic phase" [14] 
(Fermi liquid of quarks, with baryon-like excitations near the Fermi surface). Thus 
far, this conjecture was justified by qualitative arguments only, and we would like to 
make a quantitative argument by working out the details and the implication of the 
transition. 

• In the skyrmion model of large-iV c nuclear matter, Kugler and Shtrikman [19] showed 
that at low pressures, the lowest-energy configuration of is the FCC lattice of skyrmions 
with 4 different isospin orientations, while the high-pressure configuration is the simple 
cubic lattice of half-skyrmions. Later, Manton and Sutcliffe [13] described these lattices 
in terms of instantons on a T 4 torus. We would like to re-interpret their results in 
terms of instanton lattices to see where the Klebanov-Kugler-Shtrikman transition 
point (between whole-skyrmion and half-skyrmion lattices) corresponds to the popcorn 
transition from a 3D to a 4D lattice. Also, we would like to see how the chiral symmetry 
restoration in the high-density phase works in terms of the instanton lattices. 

• Finally, we would like to see if any of our ID, 2D, 3D lattices --or their generaliza- 
tions -- appear in condensed matter rather than in the holographic nuclear physics. 
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We imagine a crystalline lattice of some high-spin atoms / ions / whatever whose ori- 
entations may form all kinds of complicated patterns — perhaps non-abelian patterns 
- depending on external parameters such as temperature, pressure, or doping. Hope- 
fully, the methods developed in this article would be useful for studying such systems. 
In this case, we would end up answering questions that have yet not been asked. 
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A Energy of the Instanton Zigzag 

In this Appendix we calculate the net energy of a zigzag-shaped chain (5.2) of small instan- 
tons whose orientations follow the ansatz (5.15) from section 5, 

y n = exp(m|r 2 ) x exp(i(§ + (-l) n §)n). (5.15), 

and then we minimize that energy with respect to the orientation moduli a, /3, and <j). 

The energy of a general 2D instanton configuration in a background with M2 <C M3 ~ M4 
is spelled out in eq. (5.9): 

for 

Q g (m,n) = \ + tr 2 (^ n ) + C ?) ^ ti 2 {yly n {-iT a )) 

a=l,2 

+ (l-2C 3 )tr 2 (yiy n (-ir-N mn )), (A.2) 

where C 3 = Ml/(M 2 + Mf ), < C 3 < ±, and N^ n = (X& - X£)/\X m - X n \. For the zigzag 
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arrangement of the instanton centers, 



\X m -X n \ 2 = (n-m 2 )D 2 , 



for even n — m, (A. 3) 



but 



|iV m -X n |^ = {n-m 2 )D 2 + Ae\ 

{n - m)Dr 1 - (-l) n (2e)r 2 } for odd n - m. (A.4) 



f • N 



y/(n - m) 2 D 2 + 4e 2 



Our first task is to calculate the orientation-dependent force coefficients Q z (m,n) for the 
instantons oriented according to the ansatz (5.15). For even £ = n — m, we have 

+ / .a±(3 \ /.^ \ / .a±(3 

VrnVn = eX V{~ l ^~ Tl ) X 6Xp ^- T 2 J X exp^+J^- n 

£0 £(p £(p 

= cos — + sin — cos(« ± p) x (2T2) + sin — sin(o; ± p) x (2x3) (A. 5) 

Zi Zi Zi 

where ± = (—1)™ = (— l) m , and consequently 

Q z (m,n) = \ + 4 cos 2 — + (1 - C 3 ) x 4 sin 2 — cos 2 (a ± /3). (A.6) 

Averaging between even and odd n for fixed distance £ = n — m, we obtain 

(cos 2 (a±/3)) = - x cos 2 (o; + /3) + - x cos 2 (a - (3) 
= cos 2 a cos 2 j3 + sin 2 a: sin 2 (3 

and hence 

(Q z {£ = n-m)) = - + 4cos 2 ^ + 4(1 - C 3 ) x (cos 2 acos 2 /3 + sin 2 asin 2 /3) x sin 2 -y 



9 
2 

(A.* 



+ 2^-1 + (1 - C 3 ) x (cos 2 acos 2 /3 + sin 2 asin 2 /3)) x (1 - cos(#)). 



Therefore, the average net energy of an instanton j^m due to 2-body repulsive forces from 
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the other instantons #n at even distances £ = n — m comes up to 



So 






5AM ^ (£D) 

even tyQ v ' 



(A.9) 



/9 



N r 



5XMD 2 



7T 2 N C 

20XMD- 



x 



y - 

\even 1^0 



TV 

C 2 : 12 



+ 2^-1 + (1 -C 3 )(cos 2 acos 2 /3 + sin 2 asin 2 /3)) x 



x 



J2 l - ™ m = ^s o (0) 



Vcven £^0 



P 



(| + (-1 + (l-C 3 )(cos 2 acos 2 /3 + sin 2 asin 2 /3)) x S o (0)J 



where 



cf. eq. (5.19a). 



So(</>) 



40(7T - 0) 



7T^ 



(A.10) 



(A.11) 



For instantons at odd distances from each other we have more complicated formulae. For 
odd £ = n — m we have 



vlnVn = eX P 



— ia ± ifi 



;) n) xexp\ii-r 2 )xexpy ,, 



-Ma ± i/3 



T\ 



(A.12) 



cos — x exp(±2/5Tij + sin — x (it 2 ) x expuarij 



cos — cos /3 ± cos — sin /3 x (iri) + sin — cos a x (ir 2 ) + sin — sin a x (^3) 



where ± = (—1)" = — 1 

tr fa(-» ? ^mn)) 

and hence 



4) m . Consequently, 

= — = x (cos — sin/3 x £D — sin — coso;x2e) (A. 13) 

y/(£D) 2 + 4e 2 V 2 2 ; 



Q z (m,n) = — + 4 cos 2 — cos 2 /3 + 6*3(4 cos 2 — sin 2 /3 + 4 sin 2 — cos 2 a 

4(1 -2C 3 ) / £6 n „ £6 \ 2 

+ .a^o ttt x cos — sin a x £D — 2 sin — cos a x 2e . 



(A.14) 



Note that this Q z (m,n) depends only on £ = n — m, so there is no need to average it over 
m. Instead, it is convenient to expand the right hand side of this formula and to re-group 
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the terms according to their I dependence, thus 



QziPj = (| - sin 2 (3 + cos 2 a) + (2 - sin 2 (3 - cos 2 a) x cos(£0) 

+ (l-2C 3 )(sin 2 /3-cos 2 a) x ^g = g (A. 15) 

£2^)2 _ 4e 2 

+ (1 - 2C 3 )(sin 2 /3 + cos 2 a) x ^ 2 - ^ 2 cos(£0) 
- 2(1 - 2C 3 ) cos a sin /3 x 2 sin(£0) . 

Therefore, the net energy per instanton due to repulsion from its odd-distance neighbors 



odd 5AM ^ £ 2 D 2 + 4e 2 

odd£ 



(A.16) 



is a linear combination of five series, 



S, 



n 2 N r . 



odd 



20XMD 2 



( (| - sin 2 /3 + cos 2 a) x E x + (2 - sin 2 /3 - cos 2 a) x E 2 \ 

+ (1 - 2C 3 ) (sin 2 /3 - cos 2 a) x E ; 
+ (1 - 2C 3 )(sin 2 /3 + cos 2 a) x E, 



V 



where 

Ei = 

5 2 = 

5 3 = 
E 4 = 
S 5 = 



4D 



tt -W 2 D 2 + 4e 2 

odd£ 



4L> 2 



COS' 



-2 

4D 2 



■ 2 

4L> 2 



L^i (12 f)2 . 

odd^ 

E 

odd I 

E 



7T ^W 2 L> 2 + 4e 2 

odd£ 



2n2 



D 2 - 4e 2 



odd£ L J 



-2 

4L> 2 



E 



2£)2 

2n2 



(£ 2 D 2 -4e 2 ) xcos- 



7T ^ \£ 2 D 2 + Ae 2 Y 

odd£ L J 



TV" 



E 

odd£ 



4£De x sin £(, 
[i 2 D 2 + Ae 2 Y 



- 2(l-2C 3 )cosasin/3 x E 5 / 

(A.17) 



tanh(7re/D) 

WD) ' 
sinh((7r - 20)e/L>) 
(ne/D) cosh(7re/D) ' 

1 



cosh 2 (7re/ J D) 



cosh(20 e/D) 20 cosh((7r - 2<f>)e/D) 

X 



cosh 2 {ne /D) 



IX 



cosh(7re/D) 



sinh(20e/ J D) 20 sinh((7r - 2<j>)e/D) 

X 



cosh (rre/D) n 



cosh(ne/D) 



(A.l* 



c/. eqs. (5.19.b-f). 

Altogether, combining the interaction energies between even-distance and odd-distance 



neighbors we obtain the net average energy of an instanton in a zigzag, 

ti 2 N c 

"-■per instanton = '-'even + ^odd = onTTFrvJ ^ J~ Wi^i P) (A. 19) 



for 



T = - + (-1 + (1 -C 3 )(cos 2 acos 2 /3 + sin 2 asin 2 /3)j x E o (0) 
+ (- + cos 2 a - sin 2 fi\ x Si(e/D) 

+ (^2-cos 2 a-sin 2 /3J x S 2 (0,e/D) ( 518 ) 

+ (1 - 2C 3 )(sin 2 /3 - cos 2 a) x E 3 (e/L>) 
+ (1 - 2C 3 )(sin 2 /3 + cos 2 a) x £ 4 (0, e/D) 
- 2(l-2C , 3 )cosasin/3x E 5 (<j>,e/D). 

Or rather, this is the net zigzag's energy (per instanton) due to 2-body interactions, but 
there is also a 1-body potential energy 

£ P ot = N c \MMl x e 2 (A.20) 



-2 
iituii ueineia ucnig, at ^i 

we arrive at 



due to instanton centers being at X 2 = ±e 7^ 0. Adding this potential energy to eq. (A. 19) 



CmLnton = ^^ * ^(0, «, /3) + N C \MM* X £ 2 (A.21) 

and hence eq. (5.17). Quod erat demonstrandum. 

•k ~k ~k 

In the second part of this Appendix we minimize the zigzag's energy with respect to the 
orientation moduli (f>,a,/3. Specifically, we are going to minimize J-((f),a,f3) in two stages: 
At first we hold fixed and minimize WRT a and /3, and then we allow to vary and seek 
the overall minimum. 

For the first stage it is convenient to change variables 

a — >■ a — cos a and (3 — > b = sin/3, — 1 < a, b < +1. (A. 22) 

Then in terms of (0, a, b) the J 7 function becomes a polynomial in a and b with 0-dependent 
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coefficients, 

F{(f>,a,b) = F((f)) + G((/>) x a 2 + H{<j>) x b 2 - 2,7(0) x ab - K(<p) x a 2 b 2 (A.23) 
where 

F(<P) = \- E o (0) + ^h + 2E 2 (0), 

G(0) = (1-C 3 )S O (0) +Ei - S 2 (0) - (1-2C 3 )E 3 + (1-2C 3 )E 4 (0), 

#(0) = (1-C 3 )E O (0) -S x - E 2 (0) + (1-2C 3 )E 3 + (1-2C 3 )E 4 (0), ( A - 24 ) 

J(0) = (1-2C 3 )E 5 (0), 

K(<f>) = 2(1-C 3 )E O (0). 

Minimizing the polynomial (A.23) WRT a and 6 works differently for <fi = or cf) = n and 
for generic 0's, so let us start with the special cases. For <fi = -k we have S = S 5 = while 
S 2 = —Ei and E 4 = — E 3 , hence 

K = J = H = while G > 0. (A.25) 

Consequently, J 7 becomes independent of b, but that's OK because the j3 angle does not 
affect the instantons' relative orientations ]fi m y n for <fi = k, cf. eqs. (A. 5) and (A. 12). In 
other words, /3 or b becomes an irrelevant variable at = 7T -- just like the geographic 
longitude becomes an irrelevant coordinate at latitudes ±90°. As to the a variable, G > 
means that the minimum is at a = 0, i. e., a = ±|. Altogether, we have 

vlVn+x = ±iT 3 ( A - 26 ) 

which means the anti-ferromagnetic phase AF of the instanton orientations. The energy of 
this phase corresponds to 

T-(AF) = F(0 = vr) = 5 + l El . (A.27) 

Similarly, for cf) = we have So = S5 = while S 2 = +Ei and S4 = +S 3 , hence 

K = J = G = while H < 0. (A.28) 

This time, J 7 becomes independent of a, but that's OK since the a angle becomes irrelevant 
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to the instanton's relative orientations y\j^j n for = 0, cf. eqs. (A. 5) and (A. 12). As to the 
b variable, H < calls for maximal b 2 = sin 2 /3, thus b = ±1 and (3 = if. Altogether we 
obtain 

yly n+1 = ±ir x (A.29) 

and hence a new anti-ferromagnetic phase AF'. However, this new anti-ferromagnetic phase 
AF' has a higher energy than the old anti-ferromagnetic phase AF — 

J-(AF') = F{(j> = 0) + H{<f> = 0) = | + |ej - 2(1 - 2C 3 )S 3 > .F(AF) (A.30) 

— so it is not going to survive the minimization of energy WRT <fi. 
Now consider generic values of 0, < < 7r. This time we have 

K > 0, J > 0, G - H > 0, (A.31) 

although the signs of G and H themselves depend on and the parameters e/D and M3/M4. 
Thanks to K > 0, the polynomial (A. 23) cannot have any minima at generic values of a and 
b. More precisely, a, b are limited to the square —1 < a, b < +1, and J 7 cannot have any 
local minima in the interior of this square (i.e. for — 1 < a, b < +1) except maybe at the 
center a = b = 0. 

To see this, consider any straight line through the center of the square. Along such a 
line the ratio a/b is fixed, so J 7 becomes a quadratic polynomial of radius = a 2 + b 2 with a 
non-positive coefficient of the (radius 2 ) 2 term. As a function of an un-bounded real variable 
such a polynomial does not have any minima at all. But since the radius 2 variable is bounded 
- between zero at the center and the maximum at the edge of the square -- J 7 may have 
a minimum (or minima) at the center or/and at the edge. Applying this argument to all 
lines through the center, we find that the minimum or minima of T may lie at the square's 
boundaries a = ±1 or b = ±1, or at the center a = b = 0, but not at any other interior 
point. 

Let us consider such possible minima in more detail: 

The minimum at a = b = — i.e. a = ±|, (3 = — leads to 

vlVn+i = cos 2 ± sin 2 X ^ ( A ' 32 ) 

and hence the abelian phase of the instanton orientations. Local stability of this 



minimum requires H > and GH > J 2 , and its energy corresponds to 

T{ ) = F(<f>). (A.33) 

Minimum at the boundary b = ±1 (but generic —1 < a < +1) --i.e. (3 = ±| and 
generic a — leads to 

vlVn+i = (-l) n cos - x (±in) + sin- x (cos a x (ir 2 ) + sin a x (zr 3 )J (A.34) 

which corresponds to the non-abelian orientation phase NA2. The specific location of 
such a minimum is 

b = ±1, a = ±^-^ (A.35) 

and its local stability at |a| < 1 requires 

J \ 2 H 



G - K > J but [^zTk) > g • (A ' 36) 

The energy of this minimum is 

^(NA2) = m + F W - g( ^ w . (A.37) 

o In principle there could be a similar minimum at the other boundary a = ±1 (but 
generic 6). Specifically, such a minimum would happen at 

a = ±1, 6 = ±^4^ ( A ' 38 ) 

and its local stability at |6| < 1 would require 

J \ 2 G 



// - K > J but — > — . (A.39) 

However, thanks to G > H these two requirements contradict each other, so this type 
of minimum does not happen. 

Finally, we may have minima at the corners of the (a,b) square, at a, b = ±1 i.e., 
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a = or 7r while (3 = ±|. Such minima lead to 

yly n+1 = (-l)»cos|x(±i7i) + sin^(±zr 2 ) (A.40) 

and hence the non-abelian NAl phase of the instanton orientations. Local stability of 
this phase requires G > J + K while its energy corresponds to 

JF(NAl) = F(<f>) + G{<j>) + H((f>) - 2J(0) - K(<f>). (A.41) 

Altogether, for each generic <fi ^ 0, tc we have one or more local energy minima with 
respect to a and (3 moduli, and we have analytic formulae (A. 33), (A. 37), (A.41) for their 
energies (or rather T functions) T\ ](</>), J-[N 712] (</>), .F[iWll](0) as functions of the (ft 
modulus. The next step is to minimize these functions with respect to 0, but unfortunately 
the zero-derivative equations 

*%*-<>• **&*-*• ^"* 

are transcendental and do not have analytic solutions. In principle, we could solve these 
equations numerically and then seek the global minimum among the solutions, but since we 
do not seek high precision we have used a simpler numeric method: Covered the < <p < tt 
interval with a 1800-point grid, evaluated the T[ ], J : [NA2], and .F[iWll] for all grid 
points -- and also the J-"[AF] for <f> = 7r, - - and looked for the lowest value among these 
3 x 1800 + 1 data points. 

The phase diagram on figure 6 (on page 67) was obtained by repeating this numerical 
minimization procedure for 100 x 100 combinations of the e/D and M^/M^ parameters. 
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